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SUMMARY 
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This  report  describes  a  calculation  of  radiative  transfer  in  nitrogen. 

The  gas  is  uniform  and  in  local  thermodynamic  equilibrium.  The  temperatures 

and  densities  used  in  the  calculations  range  from  5,000°K  to  35,000°K  and 
-5  -i 

from  10  atmospheres  to  10  atmospheres  respectively. 

Radiative  transfer  due  to  the  lines  is  computed  oy  techniques  developed 
in  Part  1  of  this  series  of  reports.  Between  150  and  250  small  clusters 
of  nonhydrogenic  lines  are  treated  individually  (the  exact  number  depends 
on  the  conditions);  each  cluster  contains,  on  the  average,  3  or  4  lines. 

The  remaining  lines  are  treated  as  hydrogenic  or  accounted  for  by  the  use 
of  extensions  of  the  continuum  absorption  coefficients.  The  lines  are 
classified  into  various  groups  and  the  contribution  of  each  group  to  the 
overall  radiation  from  the  gas  is  distinguished,  me  ^eduction  in  line 
intensity  due  to  the  background  continuum  is  allowed  for.  The  intensity 
due  to  a  combination  of  up  to  18  types  of  continuum  and  pseudo-continuum 
processes  is  computed  with  exact  allowance  for  self-absorption.  The  results 
for  the  total  intensity  and  the  relative  contributions  of  selected  groups 
of  contributors  are  presented  graphically.  In  addition,  an  equivalent 
gray  gas  absorption  coefficient  is  defined  and  the  values  presented. 

The  results  of  the  calculation  demonstrate  that,  in  most  cases,  the 
atomic  and  ionic  lines  are  the  most  important  contributors.  The  exceptions 
to  tnis  occur  at  nigh  densities  (lO  ^  atmospheres),  where  the  bound-free 
continua  dominate,  and  at  low  temperatures  (c.  5000°K),  where  molecular  bands 
are  most  important.  When  the  gas  is  highly  ionized,  the  dominant  contributors 
are  a  few  strong  lines  arising  fr->m  same-shell  transitions  which  have  low- 
lying  initial  energy  levels  and  are  situated  in  the  near  ultra  violet.  A 
further  interesting  and  unexpected  result  is  that  the  intensities  of  the 
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strongest  lines  are  often  significantly  influenced  by  Doppler  broadening 
of  the  line  profiJee. 
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1.  INTRODUCTION 

This  part  of  the  work  is  concerned  with  the  calculation  of  the  specific 
intensity  of  radiation  in  a  volume  of  nitrogen  gas  where  the  properties  arc 
uniform  and  the  gas  is  in  local  thermodynamic  equilibrium.  It  is  shown  in 
Section  4.1  that  the  results  presented  can  easily  be  rescaled  to  give  the 
flux  from  a  uniform,  plane  slab  of  nitrogen  which  has  transparent  boundaries 
and  no  incident  radiation.  The  range  of  thermodynamic  states  of  the  gas  and 
path  lengths  were  chosen  to  approximately  correspond  to  conditions  produced 
by  a  space  vehicle  reentering  the  earth's  atmosphere.  Particular  interest 
was  taken  in  the  atomic  and  ionic  line  radiation  because  the  importance  of 
lines  as  a  mechanism  for  transporting  energy  under  these  conditions  has  been 
debated. 

Relevant  previous  work  which  included  line  radiation  is  that  of 

13  14 

Vorobyov  and  Norman  or  nitrogen  and  of  Allen  on  air. 

The  calculations  presented  here  are  based  on  the  techniques  developed 
in  Parts  1  and  2.  Chapter  2  discusses  and  presents  tables  of  the  data  used 
to  represent  the  lines.  Chapter  3  discusses  the  data  used  for  the  continuum 
processes.  Finally,  Chapter  4  describes  and  presents  the  calculation  and 


its  results. 
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2.  SPECTRAL  LINE  DATA 

2.1  Introduction 

In  a  hydrogenic  gas  (i.e.  where  a  single  electron  experiences  a  Coulomb 
potential)  the  spectral  lines  hava  distributions  in  both  frequency  and  oscil¬ 
lator  strength  which  may  be  expressed  as  simple  analytic  expressions.  In  a 
non-hydrogenic  gas,  however,  the  electron  fields  interact  and,  except  close 
to  a  photoionization  edge,  the  nydrogenic  distributions  do  not  even  approxi¬ 
mately  describe  twe  line  spectrum.  Since  it  is  the  non-analytic  lines  which 
are  of  most  importance  in  radiative  transfer,  the  treatment  of  lines  has  to 
be  based  on  the  spectrum  as  obtained  from  experiment  or  predicted  by 
approximate  theoretical  n.athods. 

The  frequency  distribution  of  lines  is  determined  by  the  airangement  of 
energy  levels  of  the  atom  and  ions.  Section  2.2  is  therefore  devoted  to 
describing  the  energy  level  structure  of  nitrogen  with  the  aid  of  a  number 
of  diagrams. 

Section  2.3  discusses  the  method  used  to  obtain  oscillator  strengths 
and  derives  expressions  for  the  relative  rrultiplet  strengths  in  some  cases 
for  which  values  have  not  previously  been  tabulated. 

Section  2.4  considers  the  application  of  the  methods  developed  in 
Part  2  to  the  weak  lines  of  nitrogen. 

The  treatment  of  line  pro'iles  which  was  discussed  at  length  in  Part  i, 
Chap.  4  is  summarized  in  Section  2.5. 

Section  2.6  describes  the  calculation  of  the  data  necessary  for  computing 
the  intensity  of  a  group  of  closely  spaced  lines. 

The  detailed  line  spectrum  of  nitrogen  is  extremely  complicated  and  in 
order  to  carry  out  radiative  transfer  calculations  it  is  necessary  to  reduce 
the  number  of  lines  by  combining  them  into  groups  wherever  possible.  This 
process  has  been  carried  out  for  the  first  four  spectra  of  nitrogen.  Whenever 
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approximations  have  been  made  by  combining  lines  into  groups,  the  aim  has 
been  to  restrict  the  resulting  error  in  the  total  specific  intensity  of 
the  gas  to  a  few  percent,  both  under  optically  thin  and  self-absorbed 
conditions.  Comments  on  the  approximations  are  to  be  found  in  Section  2.7 
for  lines  of  NI,  in  Section.  2.8  for  lines  of  Nil,  in  Section  2.9  for  lines 
of  NIII  and  in  Section  2.10  for  lines  of  NIV.  Detailed  tabulations  of  the 
data  assembled  for  the  grouped  lines  is  collected  in  Tables  2.3  to  2.14  at 
the  end  of  this  chapter.  The  sources  of  this  data  are  given  in  Section  2.11. 

Lines  which  are  not  listed  individually  in  the  Tables  at  the  end  of 
this  Chapter  and  which  are  not  part  of  the  pseudo-continua  (see  Section  2.4) 
were  treated  either  by  approximate  extensions  of  the  line  series  to  which 
they  belong  (used  for  the  high  lines  of  low-lying  series,  see  Section  2.3) 
or  by  the  use  of  hydrogenic  formulae  (used  for  lines  from  highly  excited 
initial  states,  again  see  Section  2.3). 

2.2  Energy  Levels 

The  energy  level  diagiams  of  the  radiating  particles  are  of  great 
importance  in  line  radiative  transfer  because  the  energy  transferred  in 
each  transition  depends  on  the  lower  energy  level  through  the  occupation 
number  and  on  the  difference  in  energy  levels  (i.e.  the  frequency)  through 
the  oscillator  strength  and  the  Planck  function.  In  addition,  the  energy 
level  structure  determines  (for  a  given  line  width)  whether  or  not  the 
line  overlaps  its  neighbors  (see  Part  1.  Section  3  H>. 

The  energy  levels  of  hydrogen  are  distributed  in  a  very  regular  manner 
such  that  differences  in  spin  and  angular  momentum  cause  splitting  which  is 
much  less  than  the  energy  differences  due  to  the  principal  quantum  number, 
n  (i.e.  total  momentum).  For  example,  states  with  n  =  2  suffer  the 
greatest  separation  which  is  about  5  *  10“ 5  eV  compared  to  a  difference 
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of  1.9  eV  between  the  n  *  2  and  n  =  3  states.  In  contrast,  the  lower 
levels  of  «  non-hydrogenic  material  are  widely  split  by  differences  in  spin 
and  angular  momentum.  The  order  of  increasing  effectiveness  in  splitting 
levels  for  a  given  electron  configuration  is  J,  L  and  S,  t  and  the  core 
term.  This  can  be  seen  for  nitrogen  by  referring  to  Fig.  1  which  shows  the 
energy  levels  of  atomic  nitrogen  over  the  range  10  eV  to  13.5  eV  and  to 
Figs.  2,  3  and  4  which  apply  to  N  ,  N*  and  respectively  over 

wider  ranges.  These  diagrams  will  be  of  assistance  in  the  detailed  discus¬ 
sion  of  the  next  paragraph. 

Differences  in  values  of  J  produce  very  small  differences  in  energy 
(typical  splitting  is  10**3  eV  to  10" 1  eV).  The  width  of  the  resulting 
multipiet  is  so  small  that  the  Planck  function  and  Boltzmann  distribution 
function  are  essentially  constant  across  the  width.  At  low  densities,  the 
line  half-width  is  much  smaller  than  the  multipiet  splitting  and  the  lines 
are  effectively  isolated.  At  high  densities,  however,  the  effective  line 
width  can  '  e  considerably  larger  than  the  multipiet  splitting  and  the  whole 
multipiet  behaves  like  a  single  degenerate  line  (see  the  discussion  in 
Part  1,  Section  3.4).  At  high  values  of  l  (and  of  n)  the  term  splitting 
can  be  almost  as  small  as  the  multipiet  splitting  but  at  lower  values  it  is 
larger  and  then  the  resulting  differences  in  Planck  and  Boltzmann  functions 
must  be  taken  into  account. 

Differences  in  t  can  produce  enough  splitting  to  span  the  gap  between 
energies  with  different  values  of  n,  the  effect  being  greatest  for  low- 
values  of  l  (and  n).  Finally,  the  splitting  due  to  different  core  terms 
can  be  very  large;  each  core  term  corresponds  to  a  term  of  the  ion  which 
results  when  the  lumping  electron  is  removed  so  that  the  core  term  splitting 
is  of  the  order  of  the  term  splitting  of  the  residual  ion.  Unlike  the  otner 
forms  of  splitting  discussed  above,  core  term  splitting  does  not  decrease 
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with  increasing  n. 

The  values  of  energy  level  used  in  this  report  were  obtained  either 

from  experimental  values  compiled  by  Moore^-  or  from  the  approximate  pre- 

2 

dictions  of  Gilmore,  occasionally  supplemented  by  other  sources  or  approxi¬ 
mate  methods,  these  are  indicated  in  Section  2.11. 

2 . 3  Calculation  of  Oscillator  Strengths 

A  general  discussion  of  the  methods  available  for  calculating  oscillator 
strengths  has  been  given  earlier  (Part  1,  Section  3.2).  This  section  gives 
further  details  of  the  method  chosen. 

We  recall  from  Part  1,  Section  3.2  that  the  oscillator  strength  for  a 
line  is 

gJfJJ'  =  M>  gLfLL. 

=  2.45  *  10'2(hv  eV)f(-L)  J(M)  o2  (2.1) 

where , 

gT  =  2J  ♦  1  and  gT  =  ( 2S+1  )(2Ltl ) . 

j  )  is  tabulated  for  all  cases  of  interest  in  F.ef.  3.  In  most  cases, 

^  )  may  also  be  found  from  Ref.  3  supplemented  by  Ref.  4.  An  important 

exception  is  the  ultraviolet  same-shell  transition  array  s2p3  -  sp4  (where, 

i  ‘j 

for  the  purposes  of  obtaining  j  ( J'<) ,  the  principal  quantum  numbers  of  the 

sub-shells  are  irrelevant).  However,  Rohrlich  (Ref.  5,  p.  446)  shows  that 

/  A* 

the  values  of  (>«  )  are  the  same  as  those  of  the  transition  array 
p3  -  sp2 ,  which  is  included  in  Allen's  tabulation. 

Another  type  of  transition  array  which  is  not  given  by  Alien  is 
s2p3  -  sp3  •  p  (where  the  point  indicates  that  one  p-eiectron  is  not  equiva¬ 
lent  because  of  a  difference  in  principal  quantum  number).  It  is  not 
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difficult  to  demonstrate  that  the  only  parentage  splitting  of  s2p3  results 
from  the  positive  or  negative  spin  of  the  passive  s-electron,  i.e.  each  term 
(2Stl>k  0£  g2p3  j,as  two  paren*  terms  ( 2S+2 an(j  2S^  0jr  Sp3#  From 

the  paper  of  Rohrlich^  an  expression  for  J (M)  can  bo  obtained  as  follows. 
Applying  Ref.  5  Eq.  (17)  to  the  transition 

p3((2S.l)L)s(2S1.l)Lnp(2S.l)L,  .  p3((2Stl)L)g2(ls)(2Stl)L 


where  =  S  *  1  and  n  >  2,  we  get 

J  (M)  -  2x3*  (2Stl)(2L,+l)(2S1+l)(2L+l)(2L-tl)  * 

*  W2(S,  y  ,S,  |  ;  S1,0)W2(0,0,L,L;0,L)W2(L,0,L’ ,1;L,1)  (2.2) 


where  W(a,b,c,d;e,f )  is  the  Racah  coefficient.  This  latter  quantity  is 
tabulated  in  a  number  of  places,  the  most  extensive  of  which  appears  to  be 
that  of  Nikiforov.**  In  the  case  where  e  is  zero,  the  Racah  coefficient 
takes  on  the  simple  form  (Ref.  6,  Eq.  (3.2)) 


W(a,b,c,d;0,f ) 


(-1) 


b+c-f 


/(2bVl)(2ctl) 


(2.3) 


where  is  the  Kronecker  delta.  The  use  of  certain  symmetry  properties 

also  given  by  Nikiforov  makes  possible  the  application  of  Eq.  (2.3)  to  each 
of  the  Racah  coefficients  appearing  in  Eq.  (2.2),  after  which  the  expression 
reduces  to 

=  (2L'tl)(2S1+l)  (2.4) 

A  third  type  of  transition  array  for  which  no  relative  multiplet  strength 
is  tabulated  is  sp4  -  sp3  •  t,  here  the  expressions  presented  by  Armstrong 

7 

et  al.  were  used,  the  Jahn  coefficient  being  transformed  to  a  Pacah  coeffi- 

0 

cient  according  to  the  identity  given  by  Edmonds  and  values  being  taken  from 
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C 

Nikiforov’s  tabulation. 

The  values  of  o2  available  have  been  discussed  in  Part  1,  Section  3.2. 
For  the  most  part  we  have  used  the  Hartree-Fock-Slater  values  tabulated  by 

9 

Kelly.  A  comparison  between  oscillator  strengths  calculated  from  these 

10 

values  and  the  more  accurate  Hartree-Fock  values  also  calculated  by  Kelly 
for  multiplets  in  the  ultraviolet  of  NI  is  given  in  Table  2.1,  also  shown 
are  values  calculated  by  Griem^  using  the  Bates-Damgaard  method.  It  can  be 
seen  that  the  Hartree-Fock-Slater  values  are  never  more  than  a  factor  of  3 
from  the  Hartree-Fock  values  and  are  usually  much  closer.  In  contrast,  the 
Bates-Damgaard  values  are  frequently  low  by  an  order  of  magnitude  and  no 
values  are  predicted  for  the  strong  lines  corresponding  to  transitions 
within  the  same  shell.  Table  2.2  compares  a  few  B-D  and  H-F-S  values  in 
the  visible  and  infrared  (no  H-F  values  are  available).  Here  the  jumping 
electron  is  in  a  more  nearly  Coulomb  field  and  both  methods  give  similar 
values. 

In  the  radiative  transfer  calculations  presented  here  Hartree-Fock 
values  of  oscillator  strength  were  used  in  the  few  cases  where  they  were 
available  otherwise  lines  from  low-lying  states  were  evaluaved  by  means  of 
Hartree-Fock-Slater  values  and  lines  from  highly  excited  states  by  the 
hydrogenic  expressions  described  in  the  next  paragraph. 

Where  the  end  states  of  a  transition  are  highly  excited,  ^he  potential 
will  be  nearly  that  of  a  Coulomb  field  so  that  the  energy  levels  and  oscil¬ 
lator  strengths  will  be  approximately  given  by  hydrogenic  expressions.  In 
addition,  the  total  radiation  from  all  such  lines  is  usually  quite  small  and 
any  errors  caused  by  using  hydrogenic  results  for  these  lines  will  only  have 


a  small  effect  on  the  overall  calculation  of  intensity.  We  therefore  consi- 

.  .  ...  r<2Sc+1>  i  ( 2G+.1 )  A2Sctn  >  ,  ( 2S+1 )  , 

der  the  transition  yl  L  Jmt  L  -  y(  i.  Jnt'  L  where  y 
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stands  for  the  electronic  configuration  of  the  core.  The  hydrogenic 
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approximations  to  the  initial  energy  level,  frequency  and  oscillator  strength 
are  respectively. 
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where  z  is  the  number  of  charges  acting  on  the  radiating  electron,  Ry  is 
the  energy  of  1  Rydberg  unit  (=  13.605  eV),  E^  is  the  energy  of  the  ioniza¬ 
tion  threshold  above  the  species  ground  state,  g  =  8core2m2  *s  the  statisti¬ 

cal  weight  of  the  lower  level,  g  s(2Sr+l)(2L_+l)  is  the  statistical 

core 

weight  of  the  residual  ion  and  the  numerical  value  of  26/3/3x  is  3.92. 

It  can  be  seen  from  Eq.  (2.7)  that  as  n  -*■  •  gf  becomes  proportional 
to  n-3.  Since  the  upper  states  cf  the  high  lines  of  any  series  experience 
a  Coulomb  field,  the  n-3  dependence  is  roughly  true  for  all  high  lines. 

This  provides  us  with  an  approximate  but  convenient  method  of  dealing  with 
the  high  lines  in  low  series  for  which  Kelly  has  not  calculated  any  values 

of  o’  (in  most  cases  these  are  lines  with  n  >  10).  We  treat  all  such 

lines  as  occurring  a^  the  ionization  frequency  and  having  weighted  oscillator 
strengths  given  by 

gf  =  constant/n^  (2.8) 


where  the  constant  is  evaluated  from  the  last  line  for  which  tiiere  is  a 
Hartrec-Fock-Siater  value  of  a2.  It  should  be  observed  that  the  error 
involved  in  using  Eq.  (2.8)  is  completely  negligible  since  thi^  app-r::ima- 
ticn  is  only  used  for  the  small  number  of  lines  for  which  no  Hartree  tock- 
Slater  value  exists  and  whlco  are  self-absorbed  but  isolated. 


1 
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2.4  Treatment  of  Weak  Lines 

In  Part  2,  Section  3.5  it  was  shown  that  all  transitions  above  a  certain 
initial  level,  m,  in  a  hydrogenic  ion  of  charge  ztl  can  be  accounted  for 
by  the  expression 

K  =  K  ~  [expf^RyA  m  )2kT)  -  1]  (2.9) 

v  v,ff 

where  is  the  absorption  coefficient  due  to  free-free  transitions  (see 

Section  3.3).  The  quantity  m  is  chosen  to  be  the  initial  quantum  number 
of  the  first  high  series  which  is  either  completely  merged  according  to 
effective  widths  or  entirely  optically  thin.  Since  we  are  using  a  hydrogenic 
representation  of  the  upper  series  and  a  hydrogenic  free-free  cross  section, 
we  can  adopt  Eq.  (2.9)  for  the  treatment  of  all  high  series. 

For  hydrogen,  it  has  also  been  shown  (Part  2,  Section  3.3)  that  high 
lines  in  low  series  can  be  treated  as  extensions  of  the  bound-free  continuum 
provided  that  they  are  either  optically  thick  and  effectively  merged,  or 
optically  thin.  The  same  analysis  performed  on  the  assumed  high-line  oscil¬ 
lator  strength,  Eq.  (2.8),  would  produce  the  hydrogenic  v-3  spectral  var¬ 
iation  but  with  a  different  general  level  cetermined  by  the  constant.  In 
most  cases,  this  would  not  match  the  (non-hydrogenic )  bound-free  cross 
section  at  the  threshold  in  either  magnitude  or  slope  because  of  the  approxi¬ 
mate  nature  of  Eq .  (2.8).  While  preserving  Eq.  (2.8)  for  any  isolated,  self- 
absorbed  lines  it  is  probably  more  accurate  to  treat  other  high  lines  by 
extending  the  non-hydrogenic  bound-free  cro'-s  section  backwards.  The  slope 
of  the  bound-free  cross  section  at  the  ionization  threshold  is  only  very 
approx  lira  tel  v  known  from  our  treatment  of  the  continuum  (see  Chap.  3)  but 
the  experience  of  calculated  iion-hydrcgcnic  bound-free  cross  sections  is 
that  the  slope  at  the  threshold  is  usually  small.  We  therefore  make  the 
extended  continuum  constant  at  its  threshold  value.  In  the  case  of  a  true 
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hydrogenic  v"3  behavior  at  the  m  =  1  threshold,  the  use  of  a  constant 
cross  section  would  be  responsible  for  a  12%  underestimate  cf  absorption 
coefficient  at  a  frequency  corresponding  to  n  =  5  and  for  a  3%  underesti¬ 
mate  at  n  *  10;  self -absorption  would  reduce  these  small  errors  even 
further . 


2.5  Line  Profiles 

The  difficult  problem  of  line  broadening  has  been  discussed  at  length 
in  Part  1,  Chap.  4.  It  was  shown  there  that  lines  with  relatively  low-lying 
upper  states  can  be  treated  as  having  dispersion  profiles  with  seni-half- 
widths  given  by  electron  impact  theory  with  occasional  additional  contribu¬ 
tions  from  resonance  broadening  (in  the  c^se  of  lines  of  NI )  and  from 
Doppler  broadening.  Part  1,  Section  4.6  considers  the  case  of  lines  with 
moderate  or  high  energy  upper  states  and  concludes  that  they  may  be  approxi¬ 
mately  treated  as  quas i-statically  broadened  (with  the  wing  prof'le 
2 .  S 

L(v  )  *v*  b/[v-VQ]  '  ).  Electron  impact  broadening  and  quasi-static  broadening 
are  both  aspects  of  broadening  by  charged  particles  which  is  called  in  this 
report  by  the  commonly  accepted  name  of  Stark  broadening. 

Ter  ease  of  reference,  the  expressions  for  line  semi-half-width  first 
given  in  1,  Chap.  4  are  repeated  again  here.  The  semi-half-width  due 

to  electron  impact  broadening  is 

w(eV)  =  A  S  /N.fkT(eV)]1/6  (2.10) 

where  A  is  either  determined  from  Crier's  results  *  (see  Part  1,  Section  u.s 
or  calculated  from  an  adiabatic  approximation  (Part  1,  Eqs.  (^.20)  ani  (“.?!)) 

The  semi-half-width  due  to  resonance  broadening  is  given  by 


v(eV)  -  [(A  )  exr?-E,/k7)  ♦  (A  )  exp(-E  /kT>]  — 

res  \  *  1  res  2  2  h. 


('V 
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where 


1/2 


(A 


res 


>  =  0.175  (-i)  gjf. 


iu  fv  .  (eV) 

1U 


(2.12) 


i  being  the  lower  state  of  the  perturbing  transition.  The  lower  states 
which  are  considered  to  be  sufficiently  low-lying  to  contribute  are 
2s2 2p3  4P,  2s2 2p3  2D  and  2s2  2p3  2P. 

The  semi-half-width  due  to  Doppler  broadening  is  given  in  Part  1, 
Section  **.2.  It  is 


vo  /2kT  0 

w  =  —  / -  in  2 

D  c  /  «R 

The  line  width  parameter  due  to  qu  si-static  broadening  is  given  by 


where 

A  =  0.691  *  It)-2  [  (n'  J2  -n2  ]z-1  (2.15) 

r  -v 

If  the  lower  state  is  not  quasi-r.taticaily  broadened,  the  value  of  n  in 
Eq.  (2.15)  is  taken  to  be  zero. 

Values  of  A  are  given  in  Tables  2.3  to  2.1^  while  the  relation  of 
these  values  to  tt.“  formulae  of  this  section  is  given  in  Section  2.11. 

2 . 6  Hultipict  Growth  Parameters 

Pai :  1,  Sections  3 . u  and  3.5  discussed  the  curve  of  growth  of  a  closely- 
spaced  group  of  lines.  Some  comments  will  be  made  here  on  three  of  the 
parameters  introduced  bv  that  discussion. 
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The  combined  cluster  width,  D,  was  obtained  by  inspection  of  the  energy 
levels  listed  by  Moore. ^  In  the  case  of  states  not  covered  by  Moore’s 
table,  the  width  of  a  similar  group  of  lines  was  used.  No  great  effort  at 
accuracy  was  made  since  the  treatment  of  the  transition  region  in  which  D 
plays  a  role  is  itself  only  an  approximation. 

it 

Values  of  the  effective  number  of  lines,  N  ,  were  calculated  with  some 
care  for  the  lower,  dispersion  profile  lines.  The  values  used  for  the  higher 
states  are  slightly  less  accurate  for  the  following  reason.  Values  of  N 
for  the  higher  states  were  first  computed  on  the  basis  of  dispersion  profiles. 
It  was  subsequently  decided  to  treat  these  lines  as  quasi-statically  broadened. 
However,  the  labor  of  recalculating  correct  values  of  N  using  the  revised 
profile  was  not  felt  to  be  justifiable.  It  turns  out  that  the  values  used 
are  slightly  greater  than  the  correct  values. 

it 

The  quantity  M  ,  the  number  of  isolated  clusters,  is  nearly  always 
unity  in  the  case  of  the  quasi-statically  broadened  lines  (where  it  is  not 
unity,  the  slight  error  described  in  the  previous  paragraph  occurs  in  this 

it 

case  also).  The  values  of  M  used  for  the  lines  with  dispersion  profiles 
were  calculated  correctly  but  the  choice  of  non-merging  clusters  is  somewhat 

it 

arbitrary  and  hence  the  values  of  M  are  also  somewhat  arbitrary. 

2 . 7  Lines  of  NI 

The  study  of  hydrogen  presented  in  Part  2  and  the  work  of  Biberman 
12  13 

et  al.  ’  have  shown  that  we  may  expect  a  few  lines  from  low-lying  states 
to  be  dominant  in  an  optically  thin  slar  but  the  higher  lines  to  oecome  of 
relatively  greater  importance  as  the  optical  thickness  increases.  When 
lines  are  optically  thick,  the  position  of  the  peak  of  the  Planck  function 
is  an  important  parametei  and  the  relative  importance  of  any  line  from  ar. 
excited  state  will  he  greatest  when  it  is  self-absorbed  and  near  the  maximum 
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of  the  Planck  function. 

We  found  for  hydrogen  that  of  two  corresponding,  self-absorbed  lines 
the  one  with  the  higher  initial  energy  level  can  sometimes  transport  more 
radiant  energy  because  the  higher  line  is  at  a  lower  frequency  ana  (at 
sufficiently  low  temperatures)  lies  nearer  to  the  Planck  maximum.  In  a 
non-hydrogenic  gas  the  same  effect  occurs  but,  in  addition,  parentage 
splitting  gives  rise  to  lines  which  have  high  energy  levels  (and  hence  low 
occupation  numbers)  but  without  any  compensating  reduction  in  frequency. 
Because  of  this  effect,  the  lines  with  the  lowest  energy  core 
(ls22s22p2  3P  for  NI)  are  the  most  important.  A  schematic  representation  of 
the  frequency  distribution  of  the  lines  of  NI,  Nil  and  NIII  is  drawn  on 
Pig .  5 , 

The  expressions  used  for  line  widths  have  been  given  in  Section  2.6. 

The  various  broadening  mechanisms  were  applied  as  follows  to  the  lines  of 
NI.  Resonance  broadening  was  allowed  for  in  the  cases  of  lines  with 
2s2p4,  2s22p2(3P)3s  or  2s22p2(3P)3d  upper  states,  details  may  be  found 
in  Tables  2.3b,  2.4b  and  2.5b.  Doppler  broadening  was  allowed  for  on  all 
2s-2p  transitions.  Stark  broadening  was  allowed  for  in  all  lines. 

It  remains  for  us  to  state  at  which  level  the  change  from  dispersion  to 
quasi-static  profiles  is  considered  to  occur.  In  this  context,  it  is  unfor¬ 
tunate  for  us  that  Griem's  calculations11  cover  fewer  lines  of  NI  than  of 
Nil  and  as  a  consequence  we  have  no  direct  evidence  of  how  to  treat  the 
states  4d  and  4f.  Since  the  energy  levels  of  NI  are  generally  more  closely 
spaced  than  those  of  Nil  we  can  expect  to  find  the  onset  of  hydrogenic 
behavior  at  a  lower  quantum  number  in  NI  than  in  Nil.  Therefore,  for 
want  of  a  better  method,  we  treat  the  4d  and  4f  states  of  NI  as  quasi- 
statically  broadened. 
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A  brief  discussion  of  the  lines  of  HI  according  to  initial  states  now 
follows  and  full  details  of  the  data  values  used  may  be  found  in  Tables  2.3 
to  2.9  at  the  end  of  this  chapter. 

Initial  Term  ls22s22p3  4S 

Lines  from  this  term  lie  in  the  vacuum  ultraviolet  and  he  ce  are  of 
greatest  importance  at  high  temperatures.  Even  at  the  highest  temperature 
within  our  range  (kT  =  3  eV),  the  Planck  maximum  lies  to  the  left  of  these 
lines  (see  Fig.  5)  and  hence,  the  most  important  lines  tend  to  be  those  with 
the  lowest  frequencies.  For  this  reason,  transitions  of  the  type  2p-nl  are 
stronger  than  those  of  the  type  2s-np.  However,  estimates  indicate  that  the 
latter  can  carry  up  to  20%  as  much  energy  as  the  2p-nl  transitions. 

For  2p-nl  transitions,  the  pseudo- continuum  is  restricted  to  principal 
quantum  numbers  higher  than  5.  All  available  Hartree-Fock-Slater  values  of 
o2  have  been  used. 

The  2s-np  transitions  are  relatively  weak  and  can  be  treated  less 
carefully.  Therefore  Hartree-Fock-Slater  values  need  be  used  for  the  first 
four  lines  only  and  the  pseudo-continuum  can  be  allowed  to  extend  back  to 
the  leading  line  if  necessary. 

A  complete  list  of  the  allowed  transitions  is  given  in  Table  3a. 

Table  3b  lists  data  for  all  lines  for  which  Hartree-Fock-Slater  values  of 
o2  were  used . 

Initial  Term  ls22s22p3  2D 

The  comments  made  about  the  lines  originating  on  the  term  2s22p3  4S 
also  apply  to  the  lines  considered  here.  The  treatment  is  identical  with 
the  exception  that  all  multiplets  with  n  =  3  or  4  and  the  same  electronic 
configuration  are  grouped  together.  Lines  with  a  2s22p2(*D)  core  also 
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occur;  these  are  treated  rather  less  carefully  than  tnose  of  (3P). 

Table  2.4a  gives  general  information  and  a  list  of  allowed  transitions. 
Table  2.4b  gives  data  for  all  lines  for  which  Hartree-Fock-Slater  values  of 
c2  were  used . 

Initial  Term  ls22s22p3  2-P. 

The  treatment  of  the  previous  two  sets,  of  lines  is  followed.  Lines 
with  a  *S  core  are  weak  and  may  be  largely  combined. 

Table  2.5a  gives  general  information  and  a  list  of  allowed  transitions. 
Table  2.5b  gives  data  for  all  lines  for  which  Hartree-Fock-Slater  values  of 
o2  were  used . 

Initial  Term  2s2p4  4P 

Lines  originating  on  this  state  have  relatively  high  values  of  frequency 
and  hence  are  most  important  at  high  temperaturas .  The  value  of  initial 
energy  is  high  and  hence  these  lines  tend  to  be  relatively  weak.  The  single 
same-shell  transition  (to  2p52P)  has  a  very  high  frequency  and  may  be 
neglected. 

Table  2.6a  gives  general  information  and  a  list  of  the  transitions 
considered.  Table  2.5b  gives  data  for  all  lines  for  which  Hartree-Fock- 
Slater  values  of  o2  were  used. 

2S  +1 

Initial  Configuration  ls22s22p2(  c  Lc>3i 

There  are  many  states  consistent  with  this  configuration  and  the 
associated  energy  levels  cannot  all  be  combined.  Nevertheless,  it  is  con¬ 
venient  to  consider  together  all  the  line  clusters  from  this  group  of 


initial  states. 
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The  set  of  transitions  considered  here  includes  a  large  number  in  which 
the  jumping  electron  remains  in  the  same  shell.  The  lines  associated  with 
these  transitions  tend  to  have  large  values  of  oscillator  strength  and  to  lie 
fairly  near  (although  usually  below)  the  maximum  of  the  Planck  function. 

Of  the  transitions  involving  different  shells,  by  far  thn  strongest  are 
those  with  3t-4l  jumps.  The  31-41  transitions  fall  into  roughly  two 
classes,  those  with  low  values  of  a2  (s-p  and  p-d)  and  those  with  high 
values  of  a2  (p-s,  d-p  and  d-f).  The  values  of  a2  of  these  two  groups 
differ  by  up  to  two  orders  of  magnitude  but  under  conditions  of  self-absorp¬ 
tion,  estimates  suggest  that  the  first  group  are  not  completely  negligible. 
However,  their  relative  weakness  enables  the  lines  to  be  lumped  into  ±e wer 
groups.  Another  factor  which  frequently  helps  in  the  grouping  of  lines  is 
the  decrease  of  term  splitting  with  rising  angular  orbital  quantum  number. 

Lines  with  upper  states  such  that  n  >  5  are  weak  and  the  upper  states 
can  be  combined.  Another  group  of  relatively  weak  lines  are  those  with  a 

core.  At  most,  these  lines  carry  about  half  as  much  energy  as  those  with 
a  3P  core.  Another  feature  of  lines  with  a  core  is  that,  as  far  as 

one  can  estimate,  they  have  the  same  frequency  as  similar  lines  with  a  3P 
core.  They  have  different  initial  energies,  however,  and  are  listed  here 
separately  from  the  3P  core  lines.  In  the  computation  of  intensities,  we 
take  advantage  of  the  coincidence  of  frequencies  and  combine  the  lines  which 
have  a  core  with  the  corresponding  lines  which  have  a  3P  core;  the 

difference  in  initial  energies  is  accounted  for  in  the  calculation  of  the 
occupation  number.  Lines  with  a  core  are  completely  negligible. 

Table  2.7a  lists  the  allowed  transitions.  Table  2.7b  gives  data  for  al± 
lines  for  which  Hartree-Fock-Sdater  values  of  a2  were  used;  this  table 
classifies  the  lines  into  same-shell  and  different  shell  transitions. 
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2S  +1 

Initial  Configuration  ls22s22p2(  c  Lc)4t 

We  can  afford  to  be  somewhat  crude  about  the  treatment  of  these  lines 
since,  due  to  their  low  occupation  numbers  and  low  frequencies  they  are  weak 
compared  to  the  3-n  transitions.  A  numerics 3  investigation  of  the  same- 
shell  transitions  shows  that  the  4-4  transitions  are  at  most  only  about 
10%  as  strong  as  the  3-3  transitions  for  temperatures  above  0.5  eV.  In 
particular,  the  very  low  frequencies  of  the  4d-4f  transitions  makes  them 
completely  negligible.  Kelly  does  not  give  values  of  er2  for  the  4f  to 
ng  transitions  and  we  will  neglect  them;  these  transitions  are  probably 
comparable  to  other  4-n  transitions  but  do  not  cn  their  own  give  a  signifi¬ 
cant  contribution  to  the  tota.;  radiation. 

The  allowed  transitions  are  between  the  same  types  of  terms  as  those 
originating  from  n  *  3  states,  the  only  exception  is  that  4f  -  ng  transi¬ 
tions  also  exist.  Since  we  are  not  including  this  additional  type  of  transi¬ 
tion  in  our  calculation,  we  will  not  list  the  individual  multiplets. 

Table  2.8  gives  details  of  the  data  used  for  all  lines  whose  oscillator 
strengths  were  calculated  from  Hartree-Fock-Slater  values  of  o2 . 

Initial  Configurations  with  n  >  5 

All  series  with  n  >  5  are  sufficiently  weak  and  sufficiently  highly 
excited  to  be  treated  as  hydrogenic.  We  use  the  following  factors  to 
approximately  allow  for  multiplet  splitting; 

N  =  125  for  3P  core  ,  N  =52.7  for  core 

and  D  =  0.15/n  . 

Same-shell  transitions  are  negligible  for  n  >  7  due  to  the  very  low 
frequencies  and  for  n  =  5  or  6  only  the  s-p  and  p-d  transitions  need 
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be  included.  The  relevant  data  is  given  in  Table  2.9. 

2.8  Lines  of  Nil 

Over  an  important  region  of  our  range  of  conditions,  the  majority  of 
nitrogen  particles  are  singly  ionized.  We  therefore  account  for  the  lines 
of  Nil  in  much  the  same  way  as  the  lines  of  NI.  It  is  expected  that  the 
relative  importance  of  lines  to  continuum  will  decrease  as  the  gas  becomes 
more  highly  ionized.  This  is  because  the  energy  levels  of  the  ionic  species 
are  more  widely  separated  than  those  of  the  atom  with  the  result  that  the 
occupation  numbers  of  the  excited  states  are  low  and  the  frequencies  of 
transitions  from  the  ground  state  are  high. 

Doppler  broadening  was  allowed  for  on  all  2s-2p  transitions.  States 
with  n  >  5  were  assumed  to  be  quasi-statically  broadened. 

Initial  Configuration  ls22s22p^ 

The  lines  originating  from  this  configuration  are  less  numerous  than 
those  from  the  2s2 2p3  configuration  of  NI  because  there  is  no  parentage 
splitting  when  a  2p  electron  jumps.  A  further  simplification  can  be  made 
by  neglecting  the  2s-np  transitions  on  the  grounds  that  they  have  extremely 
large  frequencies  ( >29  eV)  and  hence  occur  where  the  Planck  function  is  very 
small. 

In  calculating  the  electron  impact  line-widths  of  the  same-shell 
transitions  (by  adiabatic  theory,  see  Part  1,  Chap.  4),  only  the  perturba¬ 
tion  of  the  lower  state  was  considered  because  values  of  c2  for  the  transi¬ 
tions  2s2p3-2p4  are  not  available.  In  the  case  of  the  equivalent 
transitions  of  NI  this  procedure  is  approximately  correct. 

Table  2.10a  lists  the  allowed  transitions.  Table  2.10b  gives  detailed 
data  for  all  lines  for  which  Hartree-Fock-Slater  values  of  o2  were  used. 
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Initial  Configuration  2s22p3l 

The  allowed  transitions  from  these  states  are  the  same  as  those  for  the 
3P  core  of  NI,  except  that  the  multiplicity  is  always  less  by  unity.  The 
frequencies  of  these  transitions  are  more  widely  spaced  than  those  of  the 
corresponding  transitions  in  NI  and  thus  their  combination  involves  a 
greater  percentage  error.  However,  the  relatively  high  values  of.  initial 
energy  level  tend  to  make  these  lines  comparatively  weak. 

Table  2.11  gives  details  of  the  line  data  for  all  lines  for  which 
Hartree-Fock-Slater  values  of  o2  were  used. 

Initial  Configurations  with  n  >  4 

All  transitions  of  Nil  with  n  >  4  are  treated  as  hydrogenic  except  the 
n  *  4  s-p  and  p-d  same-shell  transitions  which  are  listed  in  Table  2.12. 

2.9  lines  of  NIII 

The  treatment  of  these  lines  is  made  easier  by  the  fact  that  they  should 
be  quite  closely  hydrogenic.  This  occurs  because,  in  the  most  important 
transitions,  the  jumping  electron  is  alone  and  outside  a  system  of  closed 
subshells. 

We  therefore  use  the  self-consistent  f itld-function  values  of  a2  to 
calculate  only  the  same-shell  transitions  and  the  2s2 2p  -  2s2 nl  transitions. 
Details  of  these  lines  may  be  found  in  Table  2.13.  Doppler  broadening  was 
allowed  for  in  all  the  lines  listed  in  Table  2.13. 

2.10  Lines  of  NIV 

There  are  a  significant  number  of  particles  of  N++*  present  for  condi¬ 
tions  of  interest  to  us  only  when  the  temperature  is  high  and  the  density 
is  low.  Even  under  these  conditions  it  seems  probable  that  the  lines  will 
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be  relatively  weak  on  account  pf  the  very  large  differences  between  energy 
values  (see  the  comment  in  Section  2.6). 

This  anticipated  weakness  is  fortunate  since  not  many  values  cf  cr2  are 
available.  We  use  the  Hartree-Fock  values  for  the  2s-2p  transitions  and 
for  the  series  2s2 -  2snp1P.  Of  higher  transitions,  only  those  with  a 
2S  core  are  considered  and  these  are  treated  as  hydrogen ic. 

Details  of  the  more  important  lines  are  given  in  Table  2.14.  Doppler 
broadening  was  allowed  for  in  all  lines  listed  in  Table  2.14. 

2.11  Sources  of  Data 

The  following  symbols  identify  the  sources  of  the  data  tabulated  in 
Tables  2.3  to  2.14. 

Energy  levels  and  frequencies 

2 

G.  Values  estimated  by  Gilmore, 
l  J 

HF  Values  calculated  by  Kelly10  by  the  Hartree-Fock  method. 

These  values  are  accurate  to  about  20%. 

*  Estimates  obtained  by  interpolation  between  values  in  the  same  line 

series, where  known, or  otherwise  by  drawing  a  smooth  curve  from  known 

2 

values  (or  the  estimates  of  Gilmore  )  to  the  hydrogenic  values  of 
the  nigh  lines. 

**  If  applied  to  an  energy  level,  rhis  symbol  indicates  an  estimate 

obtained  from  the  known  energy  level  cf  a  corresponding  term  with  a 
different  core  by  adding  to  the  known  value  the  difference  between 
the  energy  values  of  the  ion  terns  corresponding  to  the  two  core 
terms.  If  applied  to  a  value  of  frequency,  this  symbol  indicates 
that  the  value  is  that  of  a  similar  transition  with  a  different  core 
term.  This  method  gives  correct  values  near  the  ionization  limit. 

At  lower  energies,  comparison  of  estimated  values  with  known  values 
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indicates  a  distinct  improvement  over  assuming  hydrogenic  binding 
energies. 

Unmarked  values  were  obtained  from  experimental  results  listed  by 

1  2 
Moore  or,  in  a  few  cases,  by  Gilmore. 

Oscillator  Strengths 

HF  Hartrce-:  ock  values  due  to  Kelly. 

Unmarked  values  were  calculated  as  described  in  Section  2.3  (see  also 
Part  1,  Section  3.2)  from  Hartree-Fock-Slater  values  of  o2  given 
by  Kelly.9 

Line  width  constant  A  for  Stark  broadening  (electron  impact  or  quasi-static 

broadening) 

a  Values  calculated  from  the  adiabatic  approximation  to  electron 
impact  theory  (Part  1,  Eqs.  (**.20  and  (4.21)).  The  value 
corresponding  to  the  strongest  line  of  each  group  was  calculated. 

The  other  lines  will  have  the  same  Stark  profile  to  within  the 
accuracy  of  the  calculation. 

G  The  accurate  results  of  Griem^  were  matched  to  the  approximate 
expression,  Eq.  (2. 1C)  at  20,000°K. 

*  The  width  of  a  similar  electron  impact  broadened  line  with  another 
core  term  was  used. 

(unmarked)  These  values  were  obtained  from  Eq.  (2.14)  which  assumes 
a  quasi-static  wing  shape. 


Line  width  constants  for  resonance  broadening,  A 

res 

All  values  were  obtained  from  Eq.  (2.12). 


22 


Line  width  due  to  Doppler  broadening 

This  quantity  was  calculated  from  the  expressi  n  (2.13)  which  depends 
on  the  line  being  broadened  only  through  the  line  center  frequency.  It  is 
therefore  not  necessary  to  list  the  width  of  each  line  in  t Tables.  The 
particular  lines  which  were  considered  to  be  Doppler  broadened  are 
described  in  Sections  2.7  -  2.10. 

*  * 

Hultiplet  growth  parameters,  N  ,  N  and  D 

The  calculation  of  these  quantities  is  described  in  Section  2.6. 
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TABLE  2.1 

Comparison  of  Ultraviolet  Multiplet  Oscillator  Strengths  of  NI 


Subscript  H-f-S  indicates  values  of  o2  from 

ref. 

9 

Subscript  H-F 

indicates  f 

values  from 

ref. 

10 

Subscript  B-D 

indicates  f 

values  from 

ref. 

11 

Transition 

f  r- 

f  „ 

rH-F-S 

H-F 

B-D 

2s22p34S-2s2p44P 

0.454 

0.506 

2s22p23s4P 

0.183 

0.097 

0.0048 

2s22p32D-2s2p42P 

0.403 

0.437 

2s2p42D 

0.290 

0.314 

2s22p23s2P 

0.0742 

0.0432 

0.0045 

2s22p23s2D 

0.089 

0.0584 

0.0011 

2s22p2(3P)3d2D 

0.0081 

0.0030 

0.00032 

2s22p3(3P)3d2F 

0.0451 

0.020 

0.0089 

2s22p32P-2s2p42S 

0.137 

0.149 

2s2p42P 

0,364 

0.394 

2s2p42D 

0.139 

0.151 

2s22p2(3P)3s2P 

0.0634 

0.0382 

0.012 

2s22p2( ,D)3s2D 

0.0435 

0.0314 

0.C139 

2s22p2( 3P)3d2P 

0.0119 

0.0068 

0.00067 

2s22p2( 3P)3d2D 

0.0359 

0,0144 

0.0123 

Where  no  values  of  _ 

B-D 

are  given,  the 

Bates-Damgaard  method 

is  not 

applicable 


TABLE  2.2 


Comparison  of  Visible  Multiplet  Oscillator  Strengths  of  NI 
Subscript  H-F-S  indicates  a2  values  from  ref.  9 


Subscript  Q-D 

indicates  f  values  from  ref. 

11 

Transition 

f 

H-F-S 

fB-D 

2s22p23s4P-*s22p23p‘4D 

0.594 

0.516 

2s22p23p4P 

0.377 

0.400 

2s22p23p4S 

0.139 

0.116 

2s22p*3s2P-2s22p23p2D 

0.550 

0.558 

2s22p23p2P 

0.370 

0.355 

2s22p23p4S»2s22p23d4P 

0.779 

0.930 

2s22p23p2S-2s22p23d2P 

1.057 

0.872 

2s22p23p2D-2s22p23d2F 

0.642 

0.775 
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A  key  to  the  superscripts  of  Tables  2.3 


to  2.14  is  given  in  Section  2.11 
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TABLE  2.3a 

Lines  of  NI,  Initial  Term  ls22s22p345 
General  Data  and  Allowed  Transitions 

Initial  Energy  =  0 

Ions:  2s22p23P 

Ionization  Energies  eV:  14.548 

Upper  terms  to  which  transitions  are  allowed: 

ls22s2p44P  ls22s22p2(3P)ns4P  ls22s2p3(  5S)np4? 

ls22s22p2(3P)nd4P  ls22s2p3( 3S)npV 


2s2p35S  2s2p33S 

20.326  33.781 
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TABLE  2.3b 

Lines  of  NIf  Initial  Term  ls22s22p34S 
Detailed  Data 

Upper  states  which  are  resonance  broadened  are  as  follows 


Upper  states 
2s2p44P 

2s22p2(3P)3s4P 


(A  ). 
res  1 

0.019 

0.004 

0.002 


res' 2 


2s22p2(3P)3d4P  0.002  0 

Where  (A  ),  is  associated  with  the  lower  state  2s22p34S. 
res  1 


1 

it 

Upper  State 

hv(eV) 

Rf 

M 

M  D(eV)xl03 
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TABLE  2.4a 

Lines  of  NI,  Initial  Term  ls22s22p32D 
General  Data  and  Allowed  Transitions 


Energy  *  2.383eV 

Ions:  2s22p23P  2s22p21D  2s2p33D  2s2p31D 

Ionization  Energies:  12.165  14.063  23,600  30.040 

Upper  terms  to  which  transitions  are  allowed: 

2s2p42P  2s22p2(3P)ns2P  2s2p3( 3D)np2P 

2D  nd2P  2D 

2d  2f 

2F  2s2p3(1D)np2P 

2s22p2(1D)ns2D  2D 

nd2P  2F 

2d 

2  F 
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TABLE  2.4b 

Lines  of  NI ,  Initial  Term  ls22s22p32D 
Detailed  Data 

Upper  states  which  are  resonance  broadened  are  the  following 


Upper  States 

(A  ). 

res  1 

(A  ) 

res 

2s2p42P 

0.050 

0.024 

2s2p42D 

0.039 

0.010 

2s22p2(3P)3s2P 

0.012 

0.005 

3d2L 

0.003 

0 

where  (A  ) , 

res  1 

is  associated  with  the 

lower  state  2s 

(A  )_  is  associated  with  the  lower  state  2s22p 32P. 
res  2  K 


1 - 

Upper  State 

— 

hv(eV) 

gf 

it 

N 

it 

M 

mam 

A 

2s22p2( 3P)3s2P 

0.306 

2p  -  nl 
0.432HF 

transi 

2.60 

tions 

1.0 

ii 

0.0213° 

4s2P 

10.538 

0.094 

2,60 

1.0 

ii 

0.153  3 

2s22p2(3P)3d2L 

10.620 

0.232 

5.06 

1.93 

16 

n 

0.0530 

4d2L 

11.300 

0.404 

4.66 

1.72 

28 

0.110 

5 

11.613 

0.261 

6,60 

1.0 

33 

0.173  ; 

6 

11.750 

0.147 

6.60 

1.0 

19 

0.249 

7 

11.874 

0.091 

6.20 

1.0 

16 

0.339 

8 

11.936 

0.059 

6.20 

1.0 

11 

0.442 

9 

11.986 

0.041 

6.20 

1.0 

13 

0.560 

10 

12.017 

0.029 

6.20 

1.0 

11 

0.691  ! 

1 

2s22p2(lD)3s2D 

9.972 

0.584Hr 

3.06 

1.0 

1 

0.0213*  | 

4s2D 

** 

10.872 

0.117 

3.06 

1.0 

1 

* 

0.153 

2s22p2(1D)3d2L 

11.74CHF 

0.339^ 

6.96 

2,62 

16 

it 

0.0525 

4d2L 

it  it 

13.055 

0.431 

7.04 

2.69 

28 

0.110 

1  5 

it  it 

13,377 

0.310 

9.10 

1.0 

23 

0.173 

6 

it  it 

13.549 

0.175 

9.10 

1.0 

19 

0.249 

7 

it  it 

13.684 

0.107 

9.0 

1.0 

16 

0.339 

8 

itit 

13.775 

u.  ..  ■  .  , 

0.070 

—  ,  .  ■  - 

9.0 

1.0 

11 

0.442 

1 

3i 


TABLE  2.4b  (Continued) 


!  Upper  State 


hv(eV) 

Pj 

2p 

itit 

13.857 

0.049 

13.914 

0.034 

tions 


D(eV)xlO; 


0.560 

0.691 


2s  n 


tions 


2s2p‘42P 


2s2p3(3D)3p2L 


4p2L 

5p2L 

6p2L 


2s2p3(1D)3p2L 

4p2L 

5p2L 

6p2L 


19.38  HF 

4.37  HF 

2.59 

1.0 

13.95  HF 

3.14  HF 

3.06 

1.0 

21.419r’1 

0.514 

5.56 

2.92 

22.659°* 

0.207 

5.56 

2.92 

it 

23.000 

0.087 

7.84 

1.0 

it 

23.260 

0.042 

7.84 

1,0 

1 

it 

28.435 

I 

0.225 

5.56 

2.92 

it 

28,800 

0.088 

5.56 

2.92 

it 

29.200 

0.036 

7.84 

1.0 

29,589* 

0.018 

7.84 

1.0 

0.049 

0.070* 
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TABLE  2.5a 


Lines  of  NIf  Initial  Term  ls22s22p32P 
General  Data  and  Allowed  Transitions 


Initial  Energy  s  3.575eV 


Ionsi 


2s22p23P  2s22p21D  2s22p21S  2s2p33P  2s2p31P 


s2  o_2 1 


Ionization  Energies:  10,972  12,871 


15.025 


24.513  31.648 


Upper  terms  to  which  transitions  are  allowed: 


2s2p‘42S 


2s22p2( 3P)ns2P 
nd2P 


2s2p3( 3P)np2S 


2s22p2( JD)ns2D 


ndzS 


2s2p3(  JP)np2S 


2s22p2(1S)ns2S 


nd2D 
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TABLE  2.5b 

Lines  of  NI,  Initial  Term  ls22s22p32P 
Detailed  Data 


Upper  states  which  are  resonance  broadened  are  the  following 


Upper  State 

(A  ), 

res  1 

(A 

res 

2s2p‘42S 

0.0 

0.018 

2P 

0.050 

0.024 

2D 

0.039 

0.010 

2s22o2( 3P)3s2P 

0.012 

0.006 

2s22p2(3P)3d2L 

0.0006 

0.001 

where  (A  )  is  associated  with  the  lower  state  2s22p32D  and  (A  ) 
is  associafid  with  the  lower  state  2s22p32P.  res 


Upper  State 

hv(eV) 

sf 

ft 

N 

*  '  i 

M 

D(eV)*103 

A 

— 

2s22p2(3P)3s2P 

7.362 

2p  -  n 
0.229HF 

l  transi 

3.53 

tions 

1.0 

11 

0.0213r’’ 

4s2P 

9.346 

0.049 

0.127HF 

3.53 

1.0 

10 

0.153 

0.053 

2s22p2(3P)3d2L 

9.442 

5.71 

1.93 

7.8 

4d?L 

10.112 

0.200 

5.43 

1.93 

5.6 

0.110 

5 

10.422 

0.139 

7.34 

1.0 

28 

0.173  ; 

6 

10.558 

0.079 

7.35 

1.0 

23 

0.249 

7 

10.682 

0.049 

7.19 

1.0 

14 

0.339 

8 

10.744 

0,034 

7.00 

1.0 

12 

0.442 

9 

10,787 

0.023 

6.89 

1.0 

12 

0.560 

10 

i 

10.824 

0.018 

6.75 

1.0 

10 

0.691 

2s22p2( 1D)3s2D 

8.781. 

HP 

0.1R4nr 

2.59 

1.0 

1 

0.0213.’’ 

4s2D 

9,951jVft 

0.035 

0.005Hr 

2.59 

1.0 

1 

0.153 

2s22p2( 1D)3d2L 

7.95 

2.62 

7.8 

0.041  a 

4d2L 

11.990.,,, 

0.132 

8.10 

2.62 

5.6 

0.110 

5 

12.299** 

.0.091 

9.85 

1.0 

18 

0.173 

6 

12.460ftft 

0.052 

9.85 

1.0 

23 

0.249 

!  7 

12.586** 

0.032 

9.85 

1.0 

14 

0.339 

8 

12.650** 

0.021 

9.85 

1.0 

12 

0.442 

!  9 

12.703** 

0.014 

9.65 

1.0 

12 

0.560 

j  10 

12.735 

0.011 

9.85 

1.0 

10 

0.691 

2s22p2( 3  S ) 3 

HP 

11.551^ 

0.321 

4.29 

1.96 

3 

ft 

0.053 

4 

14.185** 

0.158 

4.15 

1.63 

3 

0.110 

5 

14.466** 

0.086 

4.15 

1  0 

3 

0.173 

6 

_ 

14.645 

0.048 

4.15 

1.0 

1 

0.249 

ft  c 
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TABLE  2.5b  (Continued) 


Upper  State 

hv(eV) 

>  H 

- ft 

M 

_ 

*  i 

M 

_ " _ L 

D(eV)*10 

i  a 

J.  .  A 

2s  -  np  Transitions 

2s2p“2S 

2P 

14.821 

17.529 

2.36?”r 

0.904 

1.93 

3.54 

1.0  ' 
1.0 

1 

2 

0.00071 

0.0005S 

2D 

12.095 

2.59 

1.0 

3 

0.00074 

!2s2p3(3P)3o2L 

22.790 

0.319 

7.11 

2.76 

10 

,0.035* 

23.392 

0.128 

7.11 

2.76  ' 

10 

p.175 

i  5r»2L 

23.796* 

0.027 

7.11 

1.0  j 

8 

0.173 

|  6p2l 

24.120 

0.015 

7.11 

i.o  ! 

i 

8 

0.249 

<2s2p3( 1  p)3p2L 

28.566 

0.136 

7.55 

2.751 

10 

i  * 

0.035. 

1  4P^L 

29.731 

0.055 

7.55 

2.75 

10 

b .  175 

1  5o2L 

30.461 

0.024 

7.56 

i.o  i 

8 

p.173 

,  6p2L 

31.072 

0.011 

_ 

7.55 

_ 

1.0  i 

8 

p.249 
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TABLE  2.6a 

Linas  of  NI,  Initial  Term  ls22s2p44P 
Ceneral  Data  and  Allowed  Transitions 


Initial  Energy  *  10.926eV 

lorn,:  2s2p35S 

Ionization  Energies:  9.422 


2s2p32D  2s2p32P  2s2p34S 

15.058  17.162  22.854 


Upper  terms  to  which  transitions  are  allowed: 


2s2p3(5S)ns4S 

2s2p3(3D)ns4D 

2s2p3( 3P)ns4P 

2s2p3(  3S)ns4S 

nd4D 

nd4D 

nd'*D 

nd4D 

4P 

4P 

.  4S 

The  transition  to  upper  state  2p^‘P  is  neglected. 
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TABLE  2.6b 

Lines  of  NT,  Tniti.'l  Term  ls^2s2n44T 


Detailed 

~ata 

upper  State 

■  ■ 

hv( eV ) 

pf 

N 

.  -  -  -  -  -f- 

D(eV )*103 

A 

2s2n3(5S)3 

8.04 

3.60 

7.4 

2.0  j 

10 

-S 

0.06“ 

4 

8.42 

1.25  ; 

7.3 

1.0  j 

10 

0.110 

5 

8.88 

0.71 

6.2 

1.0  ' 

8 

0.173 

6 

8.04 

0.4?  ; 

4.6 

i.O 

8 

0.249 

2s2p3(3D)3 

12.85 

i 

6.81 

40.0 

2.0  ; 

10 

0.06" 

4 

14. 1C 

2.87 

39.0 

i.o  : 

10 

0.110 

5 

14.52 

1.65  | 

35.3 

1.0 

8 

;  0.173 

6 

14.68 

0.95  1 

1 

32.0 

1.0 

8 

0.249 

2s2p3(3P)3 

15.65 

i 

6 . 66 

17.4 

2.0  ! 

10 

9.06“ 

4 

16.31 

2.16 

17.1 

1.0 

10 

0,110 

5 

16.62 

1.17 

14.6 

1.0 

8 

0.173 

6 

16.78 

0.67  | 

i 

12.2 

1.0  1 

8 

0.249 

2s2p3(3S)3 

21.34 

\ 

2.10 

8.3 

1.69 

10 

0.06“ 

4 

22.00 

1.25  i 

8.1 

•  n 

10 

0.11C 

5 

22.31 

0.76 

6.2 

1 

6 

0.173 

6 

22.47 

0.44  j 

4.6 

8 

0.249 

TABLE  2.7a 


2  S  +1 

Lines  of  NI,  Initial  Conf ipui'ation  ls?2s22p2(  r'  L  )3P 

c 

Allowed  Transitions 
5v>  Core  (n  >_  vB,  n*  >4) 


3s4r>-np4$ 

2 i4S-nd4P 

3d4P-n’f4D 

np4P 

n  *  s4P 

n'p4S 

np4D 

3p4P-nd4P 

n'p4P 

3s2P-np2S 

np2P 

nd4D 

n's4P 

n  'pLD 
3d4D-n ' f 4D 

np2D 

3p4D-nd4P 

n*f4F 

nd4D 

n'p4* 

nd4P 

n  'p4D 

n  's4P 

3d4F-n'f4D 

3p2F-nd2p 

3d4F-n'f4F 

n's2P 

n'*4C. 

3p2P~nd2P 

n  'd4D 

nd2P 

3d2P-n'f2D 

n»s2p 

n'o2S 

3p2D-nd2P 

n'p2P 
n  'n2D 

nd2D 

nd2F 

3d2D-n'f2D 

n's2P 

n  *  f2F 

n'p2P 

n'p2D 

3d2f-ii'f2D 

n’f2F 

Tl'f2G 

r.*D2D 

3s2D-np2P 

*D  Core  (n  >  3,  n'  >  4) 

3p2F-nd2S 

3d2S-n 

np2D 

nd2p 

n 

np2F 

nd2D 

3d2P-n 

n  's2D 

n 

T, 

3p2H-nd2P 

nd2D 

n 

nd2F 

3d2D-n 

n  's2D 

n' 

3p2F-r.d2D 

n 1 

nd2P 

n  1 

nd2C 

n' 

n's2D 

n' 

3d2r-n 

n 

n 

n 

ri 
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T/.BLE  ?,7a  (Continued) 


3s2S-np2P 


*S  Core  (n  3,  n*  >  4) 

3p2P-nd2D 

n's2S 


3d2r-n**2r' 

n'f2^ 

n'f2!! 

n*p2F 


3d2D-n'f2F 

n'p2D 


TA8LE  2.7b 


2S  +) 

Lines  of  NI,  Initial  Configuration  ls22s22p2(  C  Lc)32. 

Detailed  Data 


Same-shell  transitions 


Initial 

Term 

Initial 

Enerp,y(eV) 

Final 

Term 

hv(eV) 

rf 

it 

N 

* 

M 

D(eV)xl03 

A 

3P  Core 

3s4P 

10.335 

3p4S 

1.660 

1.67 

2.07 

1.0 

1C 

3. 0342G 

3p4P 

1.509 

4.52 

6.20 

1.0 

10 

c 

3.0270 

3p4D 

1.428 

7.14 

6.60 

1.0 

11 

n 

0.0248 

3s^P 

10.689 

3d2S 

0.913 

0,460 

1.94 

1.0 

11 

3.01243 

3p2P 

1.444 

2.22 

3.54 

1.0 

16 

3.C404G 

3p2D 

1.320 

3,30 

2.60 

1.0 

11 

n 

3.0345 

3d4S 

11.995 

3d4P 

1.017 

3.12 

2.07 

1.0 

12 

0.054Ga 

3p2S 

11.602 

3d2P 

1.368 

2.11 

1.94 

1.0 

5 

3.0660G 

;3P4  5 

11.844 

3d4L 

i 

1.175 

10.79 

11.8 

1.0 

19 

0.060  3 

3p  2P 

12.126 

3d2L 

0.905 

4.14 

5.39 

1.85 

8 

0.045  a 

:3p4D 

11.763 

i  3d4L 

1.227 

18.95 

13.9 

1.0 

44 

0.058  a 

3p2D 

12.009 

:  3d2L 

0.994 

7.69 

5.10 

1.96 

16 

0.053  3 

1 D  Core 

3s2D 

i  12.356 

,  3p2L 

1.402 

10.63 

7.84 

2.92 

9 

3.017  9 

1 3p2L 

13.758 

3d2L 

1.249 

28.7 

19.5 

2.92 

10 

0.0525° 

Different-shell  Transitions 


3  P  Core 


|3s4P 

10.335 

4p4L 

2.926 

0.192 

15.2 

2.76 

16 

0.18 

n  . 

3s2P 

10.699 

| 

4p2L 

2.641 

0.086 

7.58 

2.72 

11 

0.175 

a 

;3s 

i 

10.450 

5p 

3.261 

0,230 

22.9 

| 

1.0 

j 

60 

0.173 

i 

i 

i 

j 

6p 

3.571 

0.153 

22.9 

1.0 

60 

0.249 

!. 

7P 

3.737 

0.092 

22.9 

1.0 

40 

0.339 

j 

|3p4S 

11.995 

4S4P 

0.855 

0.615 

2.07 

1.0 

!  - 

0.091 

i 

a  ; 

3p4p 

,  I 

11.844 

1  4s4  P 

1.017 

2.19 

6.20 

1.0 

!  15 ! 

0.091 

a 

3p4D 

11.763  | 

4s“p 

1.097 

3.91  ; 

6.35 

1.0 

!  15 

0.091 

a 

'3p4L  1 

li.815 

4d4  L 

1.368. 

1.15 

27.5 

1.83 

!  22 

0.110 

i _ i 
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TABLE  2.7b  (Continued) 


Initial 
j  Term 

Initial 

Energy 

Final 

Term 

hv(eV) 

. 

■  - - 

n" 

X* 

D(eV)xio 

! 

3  A 

:  3p2s 

i 

11.602 

1.313 

0.097 

1.93 

— 

1.0 

11 

! 

!  0.153 

{  3p2P 

12.126 

4s2P 

0.796 

0.86 

3.54 

1.0 

11 

1  0.153 

'  3p2D 

12.009 

4s2P 

0.912 

1.65 

2.60 

1.0 

11 

0.153 

!  3P2L 

12.002 

4d2L 

1.680 

0.54 

12,2 

1.64 

25 

0.110 

!  3p 

1 

11.875 

5 

1.825 

0.796 

51.0 

1.0 

12 

0.173 

i 

6 

2.145 

0.268 

51.0 

1.0 

12 

0,249 

! 

i 

7 

2*306 

0.150 

51.0 

1.0 

11 

0.339 

3d4P 

13.012 

4p4L 

0.309 

2.85 

9.6 

1.0 

80 

O 

CD 

O 

3d4D 

13.020 

4p4L 

0.234 

3.58 

7.7 

1.0 

37 

0.177 

3d4F 

12.989 

4p4D 

0.261 

6.93 

6.5 

1.0 

14 

0.1 70( 

:  3d2L 

13.005 

4p2L 

0.198 

4.58 

12.3 

2.66 

12 

0.175' 

3d 

13.000 

4f 

0.692 

86.0 

33.3 

1.0 

62 

0.110 

5 

0.957 

13.95 

44.1 

1.0 

62 

0.173 

j 

1 

6 

1.158 

4.41 

44.1 

1.0 

62 

0.249 

1 

1 

7 

1.257 

1.62 

44.1 

1.0 

62 

0.339 

1 D  Core 

3s2D 

12.356 

4p2L 

2.641 

0.141 

7.7 

2.9 

3 

0.175* 

5p2L 

3.261 

0.128 

7.7 

1.0 

3 

0.173 

1 

6p2L 

3.571 

0.085 

7.7 

1.0 

2 

0.249 

i 

7p2L 

3.737**1 

0.051 

7.7 

1.0 

2 

0.339 

3P 

13.758 

4 

ft* 

1.600 

ft* 

9.5 

22.5 

3.69 

20 

0.16 

5 

1.825 

0,443 

22.5 

1.0 

12 

0.173 

6 

2.145 

0.149 

22,5 

1.0 

12 

0,249 

ft* 

7 

2.306 

0.084 

22.5 

1.0 

12 

0.339 

3d 

ft  ft 

15.001 

4p 

0.248 

91.5 

19*5 

1.0 

25 

0.13  * 

ft  ft 

4  f 

0.692 

48,0 

23,2 

1.0 

60 

0.110 

ft* 

5 

0.957 

7.79 

25.0 

1.0 

60 

0.173 

ft* 

6 

1.158 

2.45 

25.0 

l.o  : 

60 

0.249 

i 

- 

7 

1.257 

0.90 

25.0 

1.0  I 

i 

60 

0.339 
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TABLE  2.8 

Lines  of  NIf  Initial  Configuration 


2S  +1 

ls22s22p2(  C  L  )4* 


« 

Different-shell  transitions 


(3P)4s 

12.876 

5p 

1.009 

22.9 

1.0 

6 

0.173 

(3P)4p 

13.318 

5s 

0.320 

13.4 

21.0 

1.0 

6 

0.173 

5d 

0.667 

3.16 

39.0 

1.0 

12 

0.173 

( 3P)4d 

13*677 

5p 

0.209 

59.0 

34.0 

1.0 

31 

0.173 

(3P)4d,f 

13.687 

5f  fd 

0.310 

58.5 

52.1 

1.0 

31 

0.173 

(3P)4 

13.550 

6 

0.607 

22.3 

125 

1.0 

50 

0.249 

7 

0.707 

7.8 

L25 

1.0 

37 

0.339 

8 

0,769 

4.0 

125 

1.0 

37 

0.442 

(*D)4 

15.448 

5 

0.310 

84.5 

68.0 

1.0 

37 

0.173 

6 

_ 

0,620 
- ...  ■ 

48.9 

— 

68.0 

1.0 

31 

0.249 
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TADLE  2.10a 

Lines  of  Nil,  Initial  Configuration  ls22s22p2 
Allowed  transitions 


Lower  term  2s22p23P 


Upper  terms 

2s2p33S 

3P 

3D 


Lower  term  2s22p21D 


2s22pns3P 

2s2p2(4P)np5S 

nd3P 

np3S 

nd3D 

np5P 

np3P 

np5D 

np3D 

2s2p2( 2P)np3S 

npJS 

np3P 

np  lP 

np3D 

np^ 

Upper  terms 

2s2p31P  2s22pns1? 

*D  ndlP 

nd1  D 
nd1  F 

Lower  term  2s22p21S 


2s2p2(2D)np1P 
np1  D 
np*F 


Upper  terms 

2s2p31P  2s22pns 1 P 

nd1? 


2s2p2(2S)np1P 
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TABLE  2.10b 

Lines  of  Mil,  Initial  Conf iruration  ls22s22p2 
Detailed  Data 


Same-shell  transitions 


Initial 

Configu¬ 

ration 

Initial 
Energy (eV) 

Final 

Configu¬ 

ration 

hv( eV) 

N 

M" 

D(eV)x.l03 

t 

I 

A  | 

i 

2s22p23P 

0 

2s2p33S 

19.232 

3.03HF 

2.74 

1.0 

16 

ai 

0.060 

3P 

13.541 

1.54Rr 

2.19 

1.0 

16 

0.060a' 

3D 

11.435 

2,13HF 

4.71 

1.0 

16 

0.0603 

2s22p21D 

1.899 

2s2p31P 

18.776 

1.23HF 

1.0 

1.0 

0 

0.067a 

! 

, 

*D 

15.977 

3.26^ 

1.0 

1.0 

0 

0.0673 

2s22p21S 

4.052 

2s2p31D 

16.623 

0.818^ 

1.0 

1.0 

o 

0.057a 

Different  -shell  transitions 


j  2s22p23P 

0 

2s22p3s3P 

18.466 

0.797HF 

— 

2.69 

1.0 

23 

0.009  r' 

1 

!  1 

4s  3P 

24.371 

0.16 

2.69 

1.0 

20 

0.051  0 

, 

3d 

HF 

n 

23.297 

3.028 

7.84 

1.85 

18 

0.0085 

i 

4d 

26.085 

1.785 

7.84 

1.85 

16 

0.135  F 

i 

i 

5 

27.275 

0.720 

9.79 

1.0 

30 

0.087 

! 

i 

l 

6 

28.068 

0.325 

9.79 

1.0 

25 

0.125 

2s22p21D 

1.899 

2s22p3s 1 P 

16.597 

0.S03HF 

1.0 

1.0 

0 

0.0065^ 

4sJP 

22.632 

0.083 

1.0 

1.0 

0 

0.043  C’ 

3d 

21.473 

1.492HF 

2.07 

2.07 

0 

o.oi  r’ 

4d 

24.262 

0,915 

2.07 

2.07 

0 

0.18  n 

5 

25.514 

0.374 

2.69 

1.0 

0 

0.087 

6 

26.172 

0.169 

2.69 

1.0 

0 

0.125 

2s22p21S 

4.052 

2s22p3s1P 

14.444 

0.108^ 

1.0 

1.0 

0 

c 

0.00b5  ’ 

4slP 

20.478 

0.015 

1.0 

1.0 

0 

0.043  0 

3dJP 

20.758 

0.296HF 

1.0 

1.0 

0 

0.011 

4dJP 

22.148 

0.167 

1.0 

1.0 

0 

0.216  r’ 

5 

23.308 

0.068 

1.59 

1.0 

0 

0.087 

6 

24 . 073 

0.031 

1.59 

1.0 

0 

0.125 
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TARLE  2.11 

Lines  of  MIT,  Initial  Configuration  ls22s22p3?. 


Same-shell  transitions 


Initial 

Initial 

final 

hu(eV) 

Pf 

* 

N 

D(eV)*i03 

A  | 

States 

Enenry(eV)  States 

— 

3s  3P 

18,482 

3pft- 

2.393 

7.97 

9.74 

2.79 

53 

0 

.0078^ 

3s*P 

10,496 

3o‘L 

2.901 

3.086 

2.72 

2.72 

0 

0 

•0123r 

3p^S 

20.939 

3d3p 

2.485 

2.25 

2.74 

1.0 

10 

0 

.0088 

3P 

21.158 

3d 

2.132 

5.79 

7.29 

1.88 

21 

0 

9br 

3D 

20.663 

3d 

2.492 

11.28 

8.35 

1.91 

48 

0 

C  r'  '  5  i 

*  .  c 

•  Oxi.  r 

JS 

22.102 

3d1? 

1.469 

0.433 

1.0 

1.0 

0 

0 

lp 

20.408 

3d 

2.926 

2.61 

K-* 

• 

OO 

CD 

1.88 

0 

0 

.008  9r‘ 

21.598 

3d 

1.847 

2.77 

1.93 

1.93 

0 

0 

.0092  ' 

Different-shell  transitions 


i  3* 

18 

,479 

!  4P 

6.695 

1.03 

12.0  | 

5.11 

29 

0.05: 

i  3o3S 

20 

,939 

4s 

3p  i 

i 

3.449 

0.414 

2.74,' 

1.0 

21 

0.051 

:  3p3? 

21 

.159 

4s 

3P  ; 

1 

3.229 

1.161 

2,59; 

1.0 

21 

0.051 

3p3D 

23 

,245 

1  4s 

3p  ! 

.  i 

1.143 

0,685 

4.7  ! 

1.0 

21 

0.051 

1  3p!S 

22 

,102 

j  9S 

1  p 

,  r  | 

2.428 

0.097 

1.0  j 

1.0 

0 

0.0434 

3p 1  ° 

20 

,408 

: 4s 

p  ! 

4.122 

0.495 

1.0  I 

1.0 

0 

0.0434 

;  3pJD 

21 

,598 

4S 

1 P  j 

2.932 

C.586 

1.0  j 

1.0 

0 

0.0434 

!  3p3L 

20 

865 

i  4d 

3l  i 

5.182 

1.964 

42  3  ! 

2.25 

124 

0.135 

;  3p 1 L 

21 

,262 

■  4d 

lL  i 

5.108 

0.611 

4.0  ' 

2.72 

124 

0.183 

j  3d 

23 

270 

4p 

i 

l 

1.934 

9.9 

19.9  j 

1.0 

372 

0.C51 

. 

4  f 

! 

2.920 

52.1 

19.4  ‘ 

1.0 

186 

0.114 

!  3s 

18 

479 

'  rD 

. 

8.430 

0.52 

11.9  i 

1.0 

25 

0.087 

1 

op 

9.422 

0.28 

11.9  | 

1.0 

25 

0.125 

!  3p 

20 

955 

5 

6.199 

2.77 

32.0  ; 

1.0 

80 

0.087 

;  6 

7.091 

1.42 

32.0  j 

1.0 

80 

0.125 

1  3d 

23. 

,270 

5 

4.153 

10.72 

30.0  j 

l.C 

80 

0.087 

6 

4.810 

3.58 

30.0  ; 

1.0 

80 

0.125 
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TA9LF,  2.12 


'.ines  of  Nil,  Same-Shell  Transitions  with  n  =  4 


Initial! 

States 

i 

Initial 
Energy  ( eV ) 

Final 

States 

lw(eV ) 

pf 

it  S  TT 

N  M  D(eV)*103 

!  !  i 

i 

a  | 

I 

' - 

! 

j  j  i  1 

4s 

|  24.423 

4p 

j  0.744  j 

13.2 

12.5  1.0  f  371 

1  | 

O 

• 

o 

tr 

4p 

j 

25.167 

4d 

i 

|  0.930 

39.4 

!  i 

:  22.6  j  1.0  j  124 

i _ L. - 1 - 1 

0.1 5^ i 
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TABLE  2.13 
Lines  of  Hill 


Same-shell  transitions 


. 

Initial 

Initial 

Tinal  1 

hv(eV)  j 

rf  ! 

n* 

*  1 
M 

0(eV)xl0 

3'  A 

States 

Fnerf»y(eV) 

States! 

_ _L 

1 

f 

2s22p2P 

0 

t 

2s2p22s| 

16 . 24 

3.340nF 

1.93 

1.0 

22 

2.37^10‘4 

a 

i 

2p: 

18.08 

3.4*  W 

3.54 

1.0 

16 

2.5  *10“U 

a 

2°l 

! 

12.52  : 

1.28  HF; 

2.59 

l.* 

10 

'>.68*10~1’ 

a 

2s2p2uP 

7.09 

2p34s! 

i 

i 

16.03 

_  „cjir 
2.05j  j 

2.08 

1.0 

17 

1.72>‘10"4 

a 

2D 

12.5 

2p! 

16.02 

0.843HP 

2.59 

1.0 

1 

2 . 1 5*19~4 

a 

2d! 

i 

12.7 

2.124HH 

3,06 

1.0 

1 

3.21*10“4 

a 

2P 

18.0 

2d 

10.44 

o.%iMr 

3.53 

1.0 

13 

3 . 86*10"“ 

a 

2  D 1 

1 

7.08 
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3 ’  CONTINUUM  CROSS  SECTIONS 

3.1  Introduction 

The  processes  which  we  consider  in  the  continuum  are  bound-fne 

transitions  of  MI,  Nil,  NIII  and  NIV,  free-free  transitions  in  the  fields 

of  N+,  N++,  N+-t  +  and  N++++  and  molecular  band  radiation  from  N?  and  N^+. 

13  14 

Other  processes  which  occur  but  have  been  shown  by  previous  work  ’  to  be 
negligible  are  free-free  transitions  in  the  fields  of  N  and  N0  and  bound- 
free  transitions  of  the  N  ion.  A  paper  by  Norman^  suggests  that  the  lat¬ 
ter  effect  can  be  significant  at  certain  wavelengths  and  so,  when  more 
reliable  cross  sections  are  available,  it  will  probably  be  worthwhile  to 
check  its  contribution  to  the  overall  radiative  energy  transfer. 

As  mentioned  in  Part  1,  Chap.  5,  data  for  these  continuum  processes  is 
more  widely  available  in  the  literature  than  data  for  the  corresponding 
bound-bound  transitions  and  it  has  not  been  necessary  to  calculate  any 
continuum  cross  sections  for  this  work. 

Section  3.2  discusses  the  sources  of  data  used  for  bound-free  transitions. 
Section  3.3  presents  the  expressions  used  for  the  free-free  cross  sections 
and  Section  3.4  describes  the  approximate  treatment  of  the  molecular  bands. 

3 . 2  Bound- Free  Cross  Sections 

Sherman  and  Kulander1^  have  used  the  Burgess-Seaton  method  to  calculate 
cross  sections  of  18  atom  -►  ion  transitions  from  low-lying  states  of  N. 

These  include  all  the  important  transitions  from  the  2p3  states  and  all 
transitions  from  2s22p2(3P)3s  and  2s22p2(3P)3p  states,  the  cross  sectioi 
are  given  for  each  frequency  at  which  a  photoionization  threshold  occurs. 

In  our  calculations,  the  photoionization  edges  belonging  to  all  terms  with 
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with  a  {^P)3s  con^ipurat ion  are  combined  at  a  single,  average  frequency  and 
all  terms  with  a  ( JP)3p  confipuration  are  combined  at  another  frequency. 

This  procedure  involves  little  error  and  is  consistent  with  the  treatment  of 
the  corresponding  line  series.  In  addition,  the  configurations  (*D)3s  and 
(‘D)3p  are  assumed  to  have  the  same  cross  sections  as  (  P)3s  and  (  P)3p 
respectively.  The  bound-free  absorption  coefficients  of  both  cores  can  then 
be  calculated  for  these  configurations  by  summinp,  the  two  occupation  numbers. 

We  treat  all  other  bound-free  transitions  of  NI  by  the  hydropenic 
expression 


o 


v 


23  h3  1  rV3  J_ 

3/3  w2c  mV  v  V 

e 


v(u) 


(3.1) 


where  p  .(v)  is  the  bound-free  Haunt  factor  associated  with  the  initial 
n  £ 

17 

state  n,£.  Values  of  pnf(v)  are  tabulated  by  Karzas  and  Latter  ,  There 

18 

exists  an  approximate  expression  due  to  Menzel  and  Pekeris  and  corrected  by 
19 

Burpess  for  the  shell  averaped  Haunt  factor  which  is  pood  out  to  an 

electron  enerpy  of  about  10eV  *rom  the  photoionization  threshold,  thereafter 

it  is  poor  but  the  cross  section  itself  is  very  low  so  that  it  is  adequate 

20 

for  our  purposes.  In  the  form  piven  by  Peach  it  is 


f^N) 


•  1  ♦  a(-I]1/3  ♦  l(f)!/3  ♦  1/n2  [-2c(fr2/3-3/3  8(f)'1/3l 


(3.2) 


♦  1/n4  2/3  B(— ) 


v^-4/3 


where  a  =  0,1728,  8  »  -0.0496  and  is  the  threshold  frequency. 

For  the  3d  initial  states,  an  examination  of  the  results  of  Karzas 
17 

and  Latter  chows  that  "he  sneli-averaped  expression  is  seriously  in  error 
but  the  follow i,ip  simple  expression  is  within  about  10%  for  'vree  electron 
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energies  less  than  about  15eV 

C  0.7626  (hv-hvT)/Ry  <_  0.131 

g(v)  =  7  (3.3) 

^  ClO(hv-hvT)/Ry]“1  (hv-hvT)/Ry  >  0.131 

No  Burgess-Seaton  calculations  appear  to  be  available  for  Nil  or 

21 

higher  spectra  (exrept  for  one  transition  of  Nil  )  but  in  any  case,  the 

ionic  bound-free  radiation  occurs  where  the  Planck  function  is  small.  It 

will  therefore  be  sufficiently  accurate  to  use  Seaton’s  approximate 
22  o 

formula  for  the  2pq  states  (this  is  claimed  to  be  correct  to  20%  for 

positive  ions)  and  hydrogenic  results  for  more  highly  excited  states. 

+++  2  2 

The  lowest  energy  configuration  of  N  is  Is  2s  and  hence  is  not  cov¬ 
ered  by  Seaton's  formula.-  The  corresponding  photoionization  is  not  expected 
to  be  very  important  and  is  therefore  also  accounted  for  by  the  hydrogenic 
formula. 

It  is  worth  remarking  that  more  accurate  cross  sections  have  been 

calculated  by  Johnston  and  his  colleagues  at  Lockheed  Missiles  and  Space 
21  7 

Company  ’  these  calculations  are,  over  a  very  wide  temperature  range, 

the  best  predictions  which  the  current  state  of  knowledge  allows  (without 

going  to  Hartree-Fock  wave  functions).  Unfortunately,  the  only  forms  in 

7 

which  these  results  are  available  are  mean  absorption  coefficients  and 

23 

spectral  sorption  coefficients  (to  lOeV  photon  energy)  in  air.  The 
latter  cannot  be  reduced  to  an  absorption  cross  section  per  particle  since 
the  bound-free  transitions  of  nitrogen  and  all  its  ions  are  lumped 
together  and  hence  depend  on  the  species  concentrations  in  air. 
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3.3  Free-Free  Cross  Sections 

Although  the  free-free  radiation  is  often  not  negligible,  it  is  never 

dominant  and  it  is  therefore  not  necessary  to  use  the  relatively  elaborate 

24  • 

method  of  Peach  to  calculate  the  free-free  cross  sections.  We  therefore 

use  for  each  ionic  species  the  hydrogenic  expression 


2  6  z?N  H 

c  =  .  R  (V,kT) 

v,ff  3/3  ch  ( 2rrm  )3/2kT1/2v3  ff 


(3.4) 


where  gff(vtkT)  is  the  free-free  Haunt  factor  which,  for  temperatures  of 

interest  to  use,  is  close  to  unity  over  the  important  frequency  range  O.leV 

to  lOeV.  We  therefore  use  Eq.  (3.4)  with  unit  Haunt  factor. 

A  remark  on  the  accuracy  of  Eq,  (3.4)  is  to  be  found  in  Part  1, 

Section  5.3,  In  the  same  section  it  was  noted  that  an  effective  value  of 
2 

z  is  often  used  to  account,  approximately,  for  non-hydropenic  effects. 

14  2 

Following  Allen  we  adopt  the  value  =  l*5  for  singly  ionized  nitro- 

25 

gen,  as  suggested  by  the  experimental  results  of  Allen  et  al,  obtained  in 

2  2 

air.  For  the  higher  ions  we  use  z  __  =  z  . 

erf 

The  inclusion  of  the  high,  merged  (or  thin)  series  into  the  free-free 
absorption  coefficient  is  accomplished  (see  Section  2.4)  by  multiplying 
Eq,  (3,4)  by  the  factor 


$  =  1  +  -££££.  {exp(z2Ry/[(m  )2kT])  -  1}  (3.5) 

core  0 . 

1 

where  m  is  the  initial  quantum  number  of  the  first  completely  merged  (or 
thin)  series  (see  Part  2,  Section  3.5).  Multiplying  Eq.  (3.4)  by  $core* 
summing  over  the  ionic  species  considered  and  inserting  values  for  z  ,  we 
obtain  for  the  absorption  coefficient  due  to  free-free  and  pseudo  free-free 
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transitions 


K 


16ir2e6 


N 


v>ff  3/3  ch  <3™  v3 


(1.5  N  ♦  4N  <f>„  „„  .  . 

N+  3P+1D  N++  2P  +  9N  -M-+*1S+  16N  ++++*2 

N  NS 


(3.6) 


Numerically,  ^  can  be  written 


K  __  =  1.75  *  10 
v,  ff 


N 

2  e 


- 777  gffd.5N  ♦ 

(kT)i/2  (hv)3  ff  N  3P+*D 


(3.7) 


+  4N  ♦  9N  ♦  +  16J  A  N 

N++2P  N+++  5  N++++  V 


where  kT  and  hv  are  in  eV  and  Ng  is  the  number  density  of  species  S 
divided  by  the  Loschmidt  number. 


3,4  Molecular  Bands 

For  densities  less  than  about  1  atmosphere,  radiation  from  molecular 
bands  is  not  significant  above  about  10,000°K  due  to  the  dissociation  o* 
the  molecules.  Below  this  temperature,  both  N^  and  N2+  contribute  mole 
cular  band  radiation.  The  determination  of  molecular  band  absorption  coeff 
cients  has  been  the  obiect  of  continuinp,  work,  both  theoretical  and 

23 

experimental.  In  particular  the  work  of  Churchill  and  his  colleagues  at 

26  27 

the  Lockheed  Missiles  and  Space  Company  and  of  Keck,  Allen  and  others  * 
at  the  AVCO  Everett  Research  Laboratory  should  be  mentioned. 

The  detailed  structure  of  a  band  absorption  coefficient  is  extremely 
complex  but  fortunately  need  not  be  considered  within  the  accuracy  of  these 
calculations.  We  therefore  look  for  a  representation  which  contains  the 


rross  features  of  the  band  shapes  according  to  some  frequency-averapinp 

23 

procedure.  Churchill  et  al.  have  published  tables  of  absorption  coeffi¬ 
cients  for  4  bands  of  and  N  +  in  air.  These  are  conveniently  fre¬ 
quency  averaged  and,  once  the  concentrations  of  and  in  air  are 

it 

known,  can  be  reduced  to  temperature  and  frequency-dependent  cross  sec- 

26 

tions  per  molecule.  Allen  has  published  tables  and  praphs  of  the  spectral 

content  of  the  local  emission  coefficient  (i.e.  of  KB)  these  pive  more 

v  v 

details  than  we  require  but  cover  7  bands  of  and  (However,  only 

enouph  additional  information  is  available  to  obtain  the  absorption  coeffi¬ 
cient  of  6  of  these  bands  and  one  of  these  is  of  neplipible  strenpth.) 
Allen’s  results  can  be  used  to  obtain  a  very  sample,  approximate  representa¬ 
tion  and  for  this  reason  they  were  adopted,  the  details  are  discussed  next. 
26 

Allen  presents  the  spectral  content  of  7  band  systems  but,  due  to 

lack  of  an  f-number,  the  absorption  coefficient  of  the  ( Birpe-Kopf ield 

23 

#1)  cannot  be  obtained  from  his  report.  Churchill  presents  results  for 
this  system  but,  since  his  value  of  f-number  is  only  tentative,  it  was  not 
considered  worthwhile  to  combine  it  with  Allen’s  spectral  predictions;  in 
any  case  molecular  transitions  are  represented  in  the  same  frequency  inter¬ 
val  by  the  related  (Birpt-Hopf ield  02)  system.  It  is  easy  to  check 

that  the  very  low  value  of  f-number  belonpinp  to  the  N'?  (Lyman  Birpe- 
Hopfield)  system  implies  an  absorption  coefficient  more  than  4  orders  of 
mapnitude  smaller  than  that  of  the  ( Birpe-Hopfield )  system  which  occu¬ 

pies  the  same  frequency  interval,  the  former  rystem  may  therefore  be 
neplected.  The  five  systems  represented  in  these  calculations  are  shown 
in  Table  3,1, 

ii 

We  are  indebted  to  Dr.  P.  R.  Churchill  of  the  Lockheed  Missiies  and  r'pace 
Company  for  supplyinp  us  with  the  species  composition  on  which  their 
calculations  were  based. 
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TABLE  3.1 


Molecular  Band  Systems  Included  in  the  Calculations 


System 

(First  Positive) 
(Second  Positive) 

N  (Birpe-Hopfield  #2) 
N?+  (First  Negative) 
(Meinel) 


Transition  in 
Spectroscopic  Notation 

A3F+  -  B3* 

R 3rr  -  C3tr 

X1E+  -  b* 

X2E+  -  B2E+ 

X2'4  -  A2^ 


Approximate  photon 
energy  ran**e  (eV) 

0.3  -  2.7 

1.8  -  5.4 

5.5  -  14.1 

1.5  -  5.0 
0.3  -  2.25 


The  absorption  coefficient  implied  in  Allen's  results  is  the  "smeared 

27  28 

rotational  line"  model  of  Keck  etal.  which  has  been  shown  by  Patch  et  al. 

to  be  equivalent  to  a  ’*iust  overlapping  line"  model.  Although,  as  indicated 

by  the  name,  the  rotational  lines  are  not  distinguished,  considerable  detail 

o *  the  vibrational  structure  is  piven.  The  linear  (true)  absorption  coecfi- 

28 

cient  has  the  form 


(3.8) 

h/kT(v  -v))cm  ^ 
oo 

where  M  is  the  number  of  molecules  per  unit  volume  which  are  in  the  lower 

electronic  state  oF  the  transition,  ♦  is  a  dimensionless  function  which 

contains  the  vibrational-rotational  structure  of  the  band  system  and 
D( r )  1  2 

i— — j  and  v  are  spectroscopic  constants, 

ea  i  oo 
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ire*  he  hv 


m  c' 

e 


l^y  kT 


P 'r- 1  N  $  expf-  h/kT( 
ea  i 


v  -v)] 
oo 


»  ,„-3  (Xcm)  * 

1  R-q6  *  10  TkTTvT 


P(r) 


ea 


♦  exPr- 
nl 
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_6 

Allen's  graphical  nresentation  is  of  the  quantity  which  is  the 

spectral  variation  of  the  local  emission  coefficient.  Now,  the  general 
shape  of  the  curves  can  in  each  case  be  approximated  by  three  straight 
lines  on  Allen's  (semi-logarithmic )  plots.  This  approximation  loses  the 
detailed  vibrational-rotational  structure  but  preserves  the  overall  band- 
shape;  it  is  better  at  high  temperatures  than  at  low  temperatures  and  better 
for  some  band  systems  (e.p.  B-H  # 2 )  than  for  others  (e.g.  N^+  Meinel), 

The  curves  were  fitted  by  eye  to  obtain  the  constants  a  and  b  in  the 
following  equation  for  each  band-system  and  temperature: 


-6 


al  9 


-v 


a2  . 


'V 


a3  e 


-v 


X  <  X  <  X, 
o  —  —  1 


1  A  <  x2 


(3.9) 


The  spectral  content  of  the  absorption  coefficient  (see  Eq.  (3.8))  is 
given  by  $  *  X  ^  «  $X  ^/(hv)^.  The  constants  in  Eq,  (3.9)  are  therefore 
redefined  so  that  we  obtain  a  quantity  n  such  that 


^X-1cm  ^  *  n/(hv)5 


(3.10) 


and 


n  * 


-8. /hv 

f  a. 

l 

e 

hv 

> 

hv 

1 

-8  /hv 

o 

<  a. 

i 

e 

hv 

> 

hv 

2 

-B3/h» 

e 

1 

a_ 

hv„ 

> 

hv 

V.  3 

2 

i  hv2  (3.11) 

—  hv3 
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where  hv  is  in  eV  and  hv.  *  hc/X..  The  constants  o.  and  8.  are 
oresented  in  Table  3.2. 


TABLE  3.2 


Values  of  a. 

and  8. 

T  °K 

_  i  _  . 

0 

1 

2 

3 

<-  a 

1  I 

0 

5.57  *106 

7.3;  *10** 

11,000 

(  8 

0 

7.177 

4.04 

l  hv 

2,695 

2.695 

0.838 

0.310 

f a 

0 

1.981*106 

3.60  *104 

8,000 

{  8 

0 

6.621 

3.69 

|  l 

1  l  hv 

2.695 

2.695 

0.838 

0.310 

1 

f  a 

0 

3.98  *105 

i.lG  *10u 

5,000 

'  8 

0 

5.53 

2.57 

hv 
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0.838 

0.310 

(2+) 

“2 

'  n 

9.588*10'18 

1.011*108 

1.88  *1013 

11,000 

J  8 

i  * 

-243.1 

10.19 

50.8 

hv 

V 

5.39 

4.396 

3.351 

1.771 

i 

1 

!  '  a 

! 

1.573»10“34 

2.452*107 

1 . 709x10 1 3 

8,000 

{  8 

-417.3 

5.761 

51.1 

[  hv 

5.39 

4.460 

3.351 

1.771 

1 

1  "  a 

1  \ 

1.78^*10"56 

2.35  »107 

7.6f»9«10  17 

5,000 

F  ! 

,<’ 8 

-647.8 

7.565 

59.1 

1  i  hv 

5,39 

4.508 

3.351 
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2 

.3 

"  a 

1.2.1 3  *10  6 

1.752*1 C* 1 

1.063*10 1  ^ 

11,000 

<  8 

-106,85 

41.07 

181.3 

hv 

14.121 

12.523 

7.514 

5.5 

ft 

1 .889*10’4 

1.181-1011 

3.632*1023 

8,000 

<  9 

-157.8 

30.05 

255.4 

hv 

14.121 

12.638 

7.491 

5.5 

"  a 

5.595*103 

4 .453*10*  0 

2.212*1030 

5,000 

/  3 

-164.3 

33.02 

371.89 

hv 

14.121 

12.473 

7.451 

5.5 

(1-) 

2 

'  a 

3.662*10 

3.046*1 07 

1 .183*10 * ' 

11,000 

/  8 

-315.2 

6.020 

31.781 

hv 

5.00 

0,000 

3.115 

1.5 

a 

2.021-10"43 

2. 029*10** 

7.645*10  1  ® 

«,000 

,  3 

-4  54 . 0 

Q 

32.007 

•  ,  hv 

5.00 

4.012 

3.147 

1.5 

ft 

-3? 

1.048*10 

2.171-105 

4.863*10® 

5,000 

/  8 

-348.0 

-4.301 

19.435 

hv 

5.00 

4.038 

3.076 

1.5 
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4.  CALCULATION  OF  RADIATIVE  ENERGY  IN  NITROGEN 

4.1  Introduction 

The  previous  chapters  or  this  part  have  described  the  calculation  and 
assembly  of  data  for  the  calculation  to  be  discussed  in  this  chapter. 

The  basic  problem  solved  is  as  follows.  It  consists  of  finding  the 
specific  intensity  of  radiation  at  a  point  in  a  volume  of  nitrogen  gas  which 
is  in  local  thermodynamic  equilibrium  and  has  the  same  thermodynamic  state 
and  species  composition  at  all  points.  The  path  on  which  the  specific  inten¬ 
sity  is  defined  originates  on  o  transparent  boundary  and  there  is  no  radiation 
incident  on  this  boundary  from  outride.  Under  these  conditions,  the  equations 
of  radiative  transfer  (Part  1,  Eq.  1.1)  can  be  reduced  to  the  quadrature 


I  ( s ) 


B^Cl  -  exp(-Kvs)]dv 


(4.1) 


where  the  origin  of  coordinate  s  lies  on  the  boundary. 

The  calculation  and  results  will  be  discussed  and  interpreted  in  terms 
of  the  problem  stated  above.  It  is  important  to  observe,  however,  that  they 
are  also  relevant  to  the  radiative  flux  emerging  from  a  homogeneous  plane 
slab  of  nitrogen  gas.  This  may  be  seen  as  follows.  Approximate,  half-range 
equations  can  be  obtained  in  a  number  of  ways;  their  most  general  form  is 
(see  Appendix  I), 


dq 


- —  =  2wa  B  -  bq 
cnv  v 


(4.2a) 


dq 

T-—  =  2?ra  B  +  bq"  (4.2b) 

dn  v  Mv 

v 

where  q*  and  q^  are  the  half-range  fluxes,  dn^  =  K^(x)dx  is  the  incre¬ 


ment  in  optical  depth  and  a  and  b  are  constants.  Solving  Eq.  (4.2a)  for 
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the  frequency  integrated  flux  emerging  from  a  homogeneous  plane  slab  of 
thickness  l  with  transparent  walls  one  gets 

q+(t)  *  2w  M  B^Ll  -  exp(-Kvbt)]dv.  (4.3) 

Comparing  Eqs.  (4.3)  and  (4.1),  one  can  see  that 

q*(l)  =  2w  i  I(M).  (4.4) 

It  is  shown  in  Appendix  I  that  if  the  approximate  equations,  Eqs.  (4.2a)  and 
(4.2b)  are  to  be  correct  in  the  optically  thick  and  optically  thin  limits, 
then  a  and  b  must  take  the  values  1  and  2  respectively.  Eq.  (4.4) 
then  becomes 

q+(l)  =  irl(2t).  (4.5) 

Eq.  (4.5)  therefore  provides  a  means  of  obtaining  a  flux  from  the  results 
presented  in  this  chapter.  Alternative  values  of  a  and  b  can  be  used 
if  desired. 

Earlier  chapters  of  this  work  have  discussed  at  length  the  calculation 
of  Ky  for  Eq.  (4.1)  and  have  developed  techniques  for  dealing  with  the 
various  difficulties  which  a.ise.  Besides  computing  I(s),  we  examine  the 
contributions  of  the  various  processes  to  determine  which  are  the  most 
important  under  any  given  condition.  This  ./as  dc  le  oy  first  dividing  the 
lines  of  each  spectrum  into  three  classes:  Class  A  containing  a  relatively 
small  number  of  strong  lines,  Class  B  containing  the  isolated,  sell-absorbed 
lines  of  low-lying  series  and  Class  C  containing  the  isolate*-1,  self-absorbed 
lines  of  all  other  series  (if  any).  The  detailed  composition  of  these  clas¬ 
ses  is  given  in  Table  4,1  which  also  acts  as  the  key  to  the  symbols  o' 

Figs.  9  to  11. 
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Section  4,2  gives  an  outline  of  the  method  of  calculation  (some  further 
details  are  to  be  found  in  Appendix  II).  Section  4.3  presents  and  discusses 
the  results. 

4.2  The  Method  of  Calculation 

The  calculation  was  programmed  in  FORTRAN  IV  (E-Level)  for  the  IBM  360 
computer  at  Brown  University.  The  program  consists  of  a  main  program  and 
sixteen  subprograms.  A  short  description  of  tnese  is  given  in  Appendix  II. 

It  is  the  purpose  of  this  section  to  give  an  outline  of  the  method  of 
computation  and  the  form  of  the  output. 

For  given  temperature  and  density,  the  program  first  calculates  the 
species  composition  according  to  the  method  described  in  Part  1,  Section  2.3. 
The  species  composition  is  consistent  with  Figs.  4  to  7  of  Part  1  and  Fig.  13 
of  this  part.  The  program  then  calculates  the  radiative  transfer  for  each  of 
a  number  of  path  lengths  according  to  the  following  procedure. 

The  program  first  deals  with  the  lines  of  each  of  the  first  four  spectra. 
It  finds  the  initial  quantum  number,  in,  of  the  first  merged  (or  thin),  high 
series  (see  Part  2,  Section  3.5).  Values  of  m  are  limited  to  m  9. 

Next,  it  treats  the  class  B  lines,  finding  for  each  series  the  first  merged 
or  thin  high  line,  n,  (see  Part  2,  Section  3.3)  and  calculating  the  inten¬ 
sity  due  to  self-absorbed  or  isolated  low  lines.  The  quantity  n  is  limited 
to  a  minimum  value  which  depends  on  the  series.  The  influence  of  the  back¬ 
ground  continuum  (excluding  the  high  lines,  see  Appendix  II)  on  each  low- 

line  int  nsity  is  allowed  for  by  the  factor  exp(-K  s),  where  K 

V  V  |C 

o  o 

is  the  continuum  absorption  coefficient  at  the  line  center;  the  justification 
of  this  procedure  may  be  found  in  Section  3.1  of  Part  1.  A  computation  simi¬ 
lar  to  that  described  above  for  the  Class  B  lines  is  now  performed  on  the 
Class  C  lines  (in  this  case  n  >  9).  In  the  treatment  of  each  series  of 


62 


either  class,  the  pseudo-threshold  frequency  defined  by  the  first  thin  or 
merged  high  line  is  stored  according  to  a  parameter  which  relates  it  to  the 
corresponding  photo ionization  transition.  Finally  for  the  lines,  the  inten¬ 
sities  due  to  the  Class  A  lines  are  computed.  Here,  the  effect  of  high  lines 
is  included  in  the  background  continuum.  The  line  intensities  are  calculated 
by  the  methods  developed  in  Sections  3.4,  3.5  and  4.7  of  Part  1. 

After  calculating  line  intensities  and  finding  the  pseudo-continuum 
thresholds,  the  program  then  finds  the  contribution  of  all  the  continuum 
and  pseudo -continuum  processes.  It  does  this  by  performing  a  Simpson's  rule 
integration  with  respect  to  v  over  the  quantity  B^Cl-expC-K^  cs)]  where 

K  is  the  total  absorption  coefficient  due  to  ill  continuum  and  pseudo- 
v  fc 

continuum  processes. 

The  contril .it ions  from  extremely  low  and  extremely  high  frequency 
continua  are  found  analytically.  For  sufficiently  low  frequencies,  the  free- 
free  absorption  coefficient  will  be  large  enough  to  produce  black-body  radia¬ 
tion.  Therefore,  if  v.  is  the  largest  value  for  which  K  s  >>  1  and 

ir  v 

hv^/kT  <<  1,  then  the  intensity  1^  due  to  frequencies  in  the  range 

% 

0  1  v  1.  v£r  f°und  by  setting  exp(hv/kT)-l  =  hv/kT  in  the  Planck  func¬ 
tion  and  integrating  from  0  to 

Iir  =  1.67  x  1010kT(hvir)serg  sec^cm"2  (4.6) 

where  kT  and  hv^  are  in  eV.  For  sufficiently  large  frequencies,  the 
continuum  absorption  coefficient  (which  is  dominated  at  large  frequencies 
by  the  bound- free  coefficient)  will  fall  off  and  the  gas  will  become  opti¬ 
cally  thin.  We  choose  an  absorption  coefficient  fall-off  of  the  form 
_  2 

«  v  for  these  high  frequencies  (this  decay  law  approximately  represents 
the  behavior  of  nonhydrogenic  bound-free  transitions).  Then,  if  is 

the  start  of  this  optically  thin  region  and  we  also  require  hv^/kT  >>  1, 
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we  can  obtain  the  intensity  I  due  to  frequencies  in  the  range 
v  <  v  by  putting  exp(hv/kT)-l  =  exp(hv/kT)  in  the  Planck  function, 

using  the  optically  thin  result 


%  - 

and  setting  Ky  =  ic/vz  where  tc  is  the  constant  of  proportionality.  The 
result  is 

*uv  =  ^.02  X  0<skT(hvuv+kT)exp(-hv/kT)  erg  sec  ^cm  ^  (4.7) 

where  hv  and  kT  are  in  eV. 
uv 

The  quantities  and  vuv  having  been  found,  it  remains  to  compute 

the  Simpson's  rule  integration  from  to  v^.  However,  this  is  not  com¬ 

pletely  straightforward  since  the  interval  v^r  to  vuv  is  usually  large 
and  the  integrand  is  usually  uneven  so  that  an  extremely  large  number  of 
intervals  are  required  to  achieve  an  acceptable  accuracy.  In  order  to  avoid 
this,  the  interval  is  divided  into  a  number  of  sections,  each  section  being 
treated  only  as  accurately  as  the  size  of  its  contribution  to  the  overall 
integral  warrants.  An  overall  accuracy  of  approximately  1%  is  aimed  at. 

•  9 

Full  details  are  given  in  Appendix  II  in  the  description  of  subprogram  CONT. 

After  calculating  the  overall  continuum  intensity,  the  program  computes 
the  contribution  of  each  individual  continuum  process.  Since  it  is  the  con¬ 
tribution  in  the  presence  of  all  the  other  processes  which  is  required,  this 
is  done  in  each  case  by  integrating  the  continuum  Intensity  in  the  absence 
of  the  process  concerned  and  subtracting  this  value  from  the  total  continuum 
intensity.  Because  of  the  tolerance  on  the  integration  procedure,  each 
individual  contribution  is  accurate  to  only  about  2%  of  the  total  continuum 
intensity. 


64 


The  output  of  the  program  consists  of  the  total  intensity  in  various 
forms.  Th>  individual  contributions  of  groups  of  lines  (see  Table  4.1)  and 
of  the  continuum  and  pseudo-continuum  processes  are  all  expressed  as  frac¬ 
tions  of  the  total  intensity.  Additional  information  which  is  given  is  the 
ratio  of  the  intensity  of  each  group  of  lines  to  the  intensity  of  the  same 
group  in  the  absence  of  a  background  continuum  and  the  quantum  numbers 
corresponding  to  the  start  of  pseudo -bound -free  continuum  absorption  coeffi¬ 
cients.  A  typical  output  of  the  program  is  presented  as  Table  4.2. 

Finally,  we  discuss  the  steps  taken  to  check  out  the  program.  In  the 

first  place,  it  must  be  said  that  to  exhaustively  check  the  operation  of 

every  possible  path  in  such  a  complex  computation  is  not  feasible.  The 

following  tests  were  made,  however.  The  operation  of  every  subprogram  was 

checked  against  hand  calculations  by  a  series  of  specially  designed  test 

programs.  Diagnostic  output  instructions  were  included  to  check  most 

intermediate  steps.  Inevitably,  when  the  subprograms  were  incorporated  into 

the  main  calculation,  changes  were  found  to  be  necessary.  In  the  event  of 

extensive  alterations,  a  subprogram  was  retested  by  means  of  the  appropriate  test 

program.  The  calculated  composition  was  compared  with  the  results  of 

29 

Drellishak  et  al.  ,  as  described  in  Part  1,  Chapter  2,  and  the  agreement 
was  satisfactory.  Random  selections  of  line  intensities  calculated  by  the 
subprograms  within  the  main  calculation  were  checked  against  hand  calcula¬ 
tions.  The  overall  continuum  absorption  coefficient  and  the  overall  energy 

at  low  temperatures  were  checked  for  order  of  magnitude  agreement  with  the 

r  3  ^4 

results  for  air  of  Churchill  et  al.‘  and  Allen  respectively.  Finally,  by 
using  a  very  low  value  of  path-length,  it  was  possible  to  use  the  program  to 
find  values  of  the  Planck  mean  absorption  coefficient.  Values  calculated  at 
kT  s  leV  and  kT  *  2eV  were  compared  with  the  Planck  mean  absorption 
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coefficients  calculated  by  Armstrong  et  al.  which  are  the  most  accurate 
currently  available.  Table  4,3  shows  this  comparison.  It  can  be  seen  that 
our  results  are  within  25%  of  Armstrong's.  This  agreement  can  he  considered 
as  satisfactory. 

4.3  Presentation  and  Discussion  of  Results 

The  results  are  shown  on  Figs.  6  to  12  while  Figs.  13  and  14  show 

additional  information.  Figures  6  to  12  contain  the  following  information. 

Figures  6a,  6b  and  6c  show  the  quantity  I/B  as  a  function  of  path  length, 

s,  for  each  of  the  temperatures  5,000°K,  11,C00°K,  15,000°K,  20,000°K, 

-1  -3  -5 

30,000°K  and  35,000°K  at  the  three  densities  10  ,  10  and  10  times  atmos¬ 

pheric.  (A  temperature  of  11,000°K  had  to  be  used  rather  than  1C,000°K  because 
of  the  availability  of  molecular  band  data,  see  section  3.4.)  The  figures 
consist  of  curves  drawn  through  values  calculated  at  six  equally  spaced  (on 
a  logarithmic  scale)  path  lengths  over  the  range  0.1  cm  to  35  cm.  The  case 
of  s  =  10  cm  has  been  cross-plotted:  Fig.  7  shows  I/B  as  a  function  of 
temperature  while  Fig.  8  shows  I/B  as  a  function  of  density.  Calculations 
for  Fig.  7  were  performed  at  temperatures  of  5,000°K,  8,000°K,  11,000°K, 
15,000°K,  20, 000° K,  25,000°K,  30,000°K  and  35,000°K.  Calculations  for  Fig.  8 
were  performed  at  values  of  p/p^  which  are  powers  of  10  in  the  range  10  ^ 
to  10"1. 

Figures  9  to  11  display  the  distribution  of  radiative  intensity  over 
the  individual  contributors.  On  these  plots,  each  contributor  is  represented 
by  an  area;  the  value  of  the  intensity  due  to  any  individual  contributor  to 
a  given  value  of  the  abscissa  may  be  found  from  the  vertical  extent  of  the 
appropriate  area.  The  lines  of  the  various  spectra  are  highlighted  by 
distinctive  cross-hatching.  The  distribution  of  intensity  is  shown  on  Figs. 

9a  to  9e  as  a  function  of  path  length,  on  Fig.  10  as  a  function  of  temperature 
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and  on  rip.  11  as  a  function  of  density.  The  curves  were  drawn  through 
values  calculated  at  the  same  points  as  used  for  the  corresponding  plots  of 
i/B  (see  preceding  paragraph).  The  symbols  used  on  Figs.  9  to  11  are 
listed  in  Table  4.1.  However,  the  following  comment  is  necessary  regard¬ 
ing  the  continuum  contributors.  As  described  in  Section  4.2,  the  radiation 
due  to  each  continuum  contributor  was  calculated  as  the  difference  between 
the  continuum  radiation  due  to  all  contributors  and  the  continuum  radiation 
due  to  all  contributors  except  the  individual  process  in  question.  When 
the  gas  is  optically  thin,  this  is  the  same  as  calculating  the  radiation 
due  to  each  contributor  alone.  However,  when  self-absorption  occurs,  the 
problem  becomes  nonlinear  with  the  result  that  the  sum  of  the  individual 
continuum  contributions  is  less  than  the  to*ai  continuum  value.  The 
difference  appears  as  part  of  the  unlabclled  area  at  the  top  of  the  higher 
density  plots  of  Figs.  9  to  11.  By  virtue  of  the  method  of  calculation, 
the  values  presented  for  each  individual  contributor  represent  the  errors 
incurred  by  neglecting  that  contributor. 

The  last  set  of  figures.  Figs.  12a,  12b  and  12c  present  the  equivalent 
gray-gas  absorption  coefficient,  K.  This  quantity  is  defined  such  that 


Before  discussing  Figs.  6-12,  we  briefly  recapitulate  the  factors 
which  govern  the  intensity  of  a  radiating  process  (a  more  detailed  discus¬ 
sion  «•% v  bn  found  in  Par*  2).  The  occupation  number  of  the  initial  state 
is  determined  by  the  total  number  of  particles  of  the  atom  or  ion  which 
are  present  and  by  the  Boltzmann  factor  (which  favors  states  of  low  enerpv). 
The  cross  section  of  a  bound-free  continuum  process  is  independent  of  the 
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thermodynamic  state  but  that  of  a  line  is  not,  due  to  the  line  orofile. 
Increasing  temperature  and  density  tend  to  increase  the  line  width,  hence 
increasing  the  radiation  from  a  self-absorbed  line.  Finally,  the  frequency 
of  a  transition  is  important  since  it  determines  the  local  value  of  the 
Planck  function;  it  was  demonstrated  in  Part  2  that  for  some  lines,  this 
can  be  the  dominant  effect. 

Figures  6  to  8,  The  Frequency  Integrated  Intensity 

Figures  6  to  8  present  I/B,  the  radiant  energy  as  a  fraction  of 
balck-body  intensity.  This  quantity  is  the  sum  over  a  large  number  of 
intensities,  the  rates  of  change  of  which  can  be  different;  it  is  there¬ 
fore  not  surprising  that  curves  of  I/B  exhibit  a  variety  of  forms.  The 
complexity  of  the  calculation  of  I/B  makes  it  difficult  to  interpret  any 
but  the  pross  features  of  its* behavior;  such  an  interpretation  is  11  that 
is  attempted  here. 

As  a  function  of  path  length,  I/B  must  always  increase  although  the 
rate  of  growth  can  be  small  when  Doppler  broadened  lines  dominate  (see 
Fip.  6c)  or  where  black-body  radiation  is  approached  (see  F ip.  6o).  Black- 
body  radiation  is  achieved  only  at  the  highest  density,  o/0q  *  10*\  and 
then  only  for  relatively  hiph  temperatures  and  path  lengths.  The  intensity 
at  5,000oK  (which,  if  may  be  noted,  has  been  multiplied  by  a  factor  of  101*) 
prows  linearly  with  s  becarre  optically  thin  molecular  bands  dominate 
the  rad;ation.  (Atomic  lines  occur  cut  are  self-absorbed  and  negligible 
compared  to  the  molecular  bands.)  The  intensity  curves  for  higher 
temperatures  prow  less  rapidly  due  to  s#lr-absorption. 

Figure  7  shows  the  quantity  (I/B )/{o/cq)  as  a  function  of 
temperature.  The  factor  (o/o  )~*  eliminates  some  of  the  density 
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dependence  of  these  curves  which  would  coincide  in  the  absence  of  self¬ 
absorption  and  change  of  species.  In  general,  the  curves  of  Fig,  7  exhibit 
an  increase  of  intensity  with  temperature.  At  o/oq  *  10  ^  this  rise  is 
due  to  an  increase  in  the  bound-free  continuum  produced  by  the  increasing 
importance  of  transitions  from  excited  states.  At  lower  densities,  the 
following  explanation  for  the  increase  of  intensity  with  temperature  is 
offered.  At  o/0q  *  10*^,  the  lines  are  not  heavily  self-absorbed.  As 
the  temperature  increases,  each  line  spectrum  in  turn  becomes  important 
and  then  dies  away  (see  Fig.  10),  Now,  the  dominant  lines  in  each  of  the 
second  and  third  spectra  (namely  TTA1,  IIIA1  and  IIIA?)  have  approximately 
the  same  frequencies  and  total  oscillator  strengths  but  at  temperatures 
where  the  third  spectrum  is  dominant,  the  Planck  maximum  is  nearer  to  the 
line  positions  than  at  temperatures  for  which  the  second  spectrum  is 
dominant.  Consequently,  the  intensity  can  increase  with  temperature.  At 

o/o  *  10  \  on  the  other  hand,  the  same  effect  is  obscured  by  the  occur- 
o 

rence  of  considerable  sel*-absorption  which  favors  the  lines  of  the  lower 
ionized  species  (because  Stark  broadening  is  more  important  for  such  lines! 

_3 

It  seems  likely  that  the  local  minimum  of  the  *  1°  curve  (which  does  not 

occur  at  °,0Q  *  10  is  due  to  th»  increased  self- absorption  of  the 
higher  species  lines. 

Tigure  8  shows  the  intensity  as  a  function  of  density.  The  rise  in 
intensity  with  increasing  density  is  due  to  the  accompanying  increase  in 
optical  thickness. 

Figures  9  to  11,  The  Pistribution  of  Padiation  Among  individual  fontriNjtors 

Figures  9a-9e  show  the  effect  of  changes  in  path  length  s  *or  a 
variety  of  conditions.  In  general,  one  '•an  see  that  the  important*  of 
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lines  with  respect  to  continuum  decreases  with  increasing  s.  This  is  due 
to  one  or  both  of  the  following  two  effects.  At  low  densities,  the  conti¬ 
nuum  is  optically  thin  while  many  of  the  stronger  lines  are  self-absorbed. 
In  this  case,  the  energy  carried  by  the  continuum  grows  linearly  with  s 
but  the  line  intensities  grow  only  as  s  ,  or  more  slowly  This  effect 

is  mainly  responsible  for  the  decrease  in  relative  line  intensity  at 
-3 

p/oq  -  10  .  The  second  effect  occurs  at  higher  densities  where  the  con¬ 

tinuum  experiences  some  self-absorption.  When  this  occurs,  the  continuum 

reduces  i-he  intensity  of  each  line  by  the  factor  exp(-K^  s)  (see  Tart  1, 

o’ 

Section  3.1),  hence  reducing  the  relative  importance  of  the  lines.  As  an 

illustration  of  the  magnitude  of  this  effect,  for  the  case  T  =  20,000°K, 

-1 

p/p  =  10  and  s  =  10rmf  the  factor  exp(-K  s)  is,  on  the  average, 

O  V  |C 

o* 

0.1.  It  is  this  effect  which  is  responsible  for  the  dramatic  fall  of  the 

line  contribution  with  s  at  p/p  -  10  In  addition  to  this  change  in 

o 

the  balance  of  energy  between  lines  and  continuum,  the  distribution  over 
line  contributors  chanre3  due,  primarily,  to  a  decrease  in  importance  as 
s  increases  of  the  Doppler  broadened  lines.  For  example,  if  two  or  more 
species  are  contributing  line  radiation,  the  relative  importance  of  the 
atom  or  lowest  ion  tends  to  increase  with  s.  This  occurs  because  the 
lines  of  higher  ions  are  predominantly  Doppler  broadened  and  hence  the 
associa*?d  line  intensity  grows  more  slowly  than  that  due  to  the  disper¬ 
sion  profiles  of  the  atomic  (or  first  ion)  lines. 

figure  10  shows  the  influence  of  temperature  on  the  distribution  of 


radiation.  At  p/p  *  10  ,  only  the  first  spectrum  makes  any  substantial 

o 

contribution  to  the  radiation.  The  lines  contribute  only  15%  or  less  of 


the  total  except  between  9,000°K  and  15,000°K  where  the  line  radiation 
contribution  rises  to  s  maximum  cf  30%.  This  maximum  is  caused  by  the 
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local  importance  of  the  IA3  group  and  this  importance  follows  from  the 

relatively  large  widths  of  these  lines.  At  higher  temperatures,  they 

become  insignificant  due  to  a  rise  in  continuum  absorption  coefficient  at 

-3 

the  frequencies  of  these  lines.  The  plot  for  o/pq  =  10  shows  the 
contribution  of  many  different  groups.  After  the  initial  dominance  of  the 
molecular  bands,  first  the  lines  and  continuum  of  the  first  spectrum  and 
then  the  lines  and  continuum  of  the  second,  spectrum  become  dominant.  In 
each  case,  the  lines  rise  into  prominence  first  and  then,  at  sliphtly 
higher  temperatures,  the  bound-free  continuum  rises  with  the  result  that 
the  relative  importance  of  the  1 ines  shows  a  number  of  maxima  and  minima. 
That  the  lines  of  a  spectrum  should  become  important  at  lower  temperatures 
than  the  continuum  is  presumably  caused  by  the  movement  of  the  maximum  of 

the  Planck  function  in  the  direction  of  higher  frequencies  as  the  tempcra- 

-5 

ture  rises.  TIk  plot  for  p/o  =  10  shows  many  of  the  features  of  that 

-3 

for  p/oq  s  10  but  the  lines  are  not  so  heavily  self-absorbed  with  the 
result  that  the  continuum  is  unimportant  except  at  low  temperatures.  This 
decrease  in  self-absorption  also  results  in  the  2s-2p  same-shell  transi¬ 
tions  making  an  appearance  in  the  first  spectrum  (group  IA1)  and  becoming 
dominant  in  the  second  spectrum  (group  IIA1).  These  lines  have  large  f- 
numbers  and  moderate  frequencies  but  suffer  from  small  Stark  widths  so 
that  they  tend  to  be  heavily  self-absorbed  at  all  but  the  lowest  densities. 
A  comparison  can  be  made  between  Fig,  10  and  Fig,  13.  Figure  13  shows  the 
distribution  of  radiating  particles  (i.e.  molecules,  atoms  and  molecular 
and  atomic  ions)  in  a  manner  similar  to  that  used  in  Fig.  10  for  the 
distribution  of  radiative  energy.  The  comparison  of  these  two  figures 
shows  tbit  the  relative  importance  of  radiative  transfer  by  atoms  and  the 
lower  ions  tends  to  lag  (in  terms  of  increasing  temperature)  behind  the 
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relative  importance  in  number  densities.  In  the  case  of  higher  ions,  the 
temperature  at  which  each  ion  is  dominant  in  the  transfer  of  radiative 
enerpy  approximately  coincides  with  the  temperature  at  which  it  is  numeri¬ 
cally  dominant.  We  may  also  point  out  that  whereas  N^+  is  not  suffi¬ 
ciently  abundant  to  appear  on  Fig.  13,  it  can  be  the  dor  inant  contributor 
to  radiative  enerpy  transfer  as  shown  on  Fig.  10;  an  explanation  for  this 
behavior  is  included  in  the  next  paragraph. 

The  influence  of  varying  density  is  shown  on  Figs.  11a  to  lid. 
Increasing  density  has  two  general  effects,  it  increases  the  number  density 
of  each  species  (and  hence  the  optical  depths)  and  it  changes  the  composi¬ 
tion  of  the  gas,  making  the  less  highly  ionized  particles  more  abundant. 

In  the  case  of  lines,  both  these  effects  favor  the  lower  spectra.  (The 
increase  in  optical  depth  favors  the  lower  spectra  because, as  pointed  out 
previously,  such  lines  tend  to  have  broader  profiles  and  hence  a  more 
rapidly  rising  intensity  than  lines  of  higher  spectra.)  Figure  11  shows  a 
fall  in  the  importance  of  lines  with  density  due  to  self- absorption  and  the 
increase  in  importance  of  the  lower  spectra  which  has  just  been  discussed. 

The  5,000°K  plot  of  Fig.  11a  shows  the  interesting  result  that  the 
intensity  from  the  N^+  molecular  bands  exceed  the  intensity  from  the 
molecular  bands  at  low  densities.  This  occurs  despite  the  number  of  parti¬ 
cles  of  N^+  being  less  than  the  number  of  particles  of  N2  by  several 
orders  of  magnitude  (compare  Fig,  13),  The  reason  for  this  behavior  is  as 
follows.  Both  bands  of  N2+  originate  from  the  ground  state  and  have 
frequencies  which  put  them  at  or  near  the  maximum  of  the  Planck  function 
at  a  temperature  of  5,000°K.  On  the  other  hand,  two  bands  of  N2  originate 
from  excited  states  which  have  low  occupation  numbers  at  this  temperature. 
The  remaining  band  of  the  Birge-Hopfield  #2  band,  originates  from  the 

ground  state  but  occurs  at  high  frequencies  where  values  of  the  Planck 


72 


function  are  low. 

We  will  now  draw  attention  to  some  additional  conclusions  which  can 
be  drawn  from  Figs.  9,  10  and  11  and  make  some  ge'~ral  comments.  From 
these  figures  it  can  be  seen  that  whereas  the  contribution  from  the  lines 
of  the  first  spectrum  is  scattered  over  many  different  groups,  one  or 
occasionally  (in  the  case  of  the  second  spectrum)  two  groups  transport  all 
the  significant  line  energy  for  the  ions.  This  is  because  most  of  the 
ion  lines  have  either  very  high  energy  levels  or  very  high  frequencies 
(which  puts  them  in  the  tail  of  the  Planck  function)  and  hence  are  very 
weak.  The  few  strong  lines  of  each  spectrum  are  the  same-shell  transitions 
from  the  ground  state.  The  balance  in  importance  between  lines  and  con- 
tinua  is  maintained  for  these  spectra  (and  if  anything,  shifted  in  favor 
of  the  lines)  because  although  few  lines  are  important,  the  bound-free 
continua  occur  at  very  high  frequencies  and  are  thus  poor  carriers  of 
radiation. 

From  Figs.  9  to  11  it  can  La  seen  that  for  temperatures  above  10,000^, 

molecular  bands  are  not  significant.  At  temperatures  above  10,000°K,  the 

-4 

lines  are  dominant  for  densities  below  approximately  10  times  atmospheric, 

-4  -2 

a  balance  between  lines  and  continua  exists  for  densities  of  10  to  10 
while  at  higher  densities  the  continuum  processes  are  the  most  important, 
although  lines  are  not  negligible. 

In  the  range  of  conditions  examined,  all  the  first  four  line  spectra 
can  be  important  but  only  the  continuum  processes  of  the  first  two  spectra 
need  be  accounted  for.  Class  C  lines  are  never  significant  while  the  only 
important  continuum  and  pseudo-continuum  mechanisms  (apart  from  molecular 
bands)  are  the  bound-free  and  hiph  lines  of  the  first  and  second  spectra. 

Finally,  it  can  be  seen  that  the  continuum  is  approximately  optically 

-2 

thin  for  densities  below  10  times  atmospheric  while  at  least  some  of  the 


73 


lines  suffer  self-absorption  at  ali  conditions  (evidence  for  this  statement 
will  be  given  in  the  discussion  of  the  equivalent  gray-gas  absorption 
coefficient). 


Figure  12 t  The  Equivalent  Gray  Gas  Absorption  Coefficient 

When  the  gas  is  optically  thin  with  respect  to  every  contributor,  K, 
defined  by  (see  Eq.  (4,8)) 

K  =  -  fcn(l  -  I/B)/s, 


is  equal  to  K^,  the  Planck  mean  absorption  coefficient.  As  I/B  -*■  1, 
small  errors  in  the  calculation  of  I  cause  K  to  become  indeterminate. 
Because  of  this  indeterminacy,  the  curves  of  K  on  Fig.  12a  for 
T  2.  20 ,000°K  could  not  be  drawn  for  high  values  of  s  and  the  curve  cor¬ 
responding  to  15 ,000°K  is  not  reliable  in  the  dashed  region. 

It  can  be  seen  from  Fig,  12  that  i<  decreases  with  increasing  path 
length  except,  possibly,  where  the  radiation  approaches  that  of  a  black- 
body  and  at  5,000°K.  In  the  5,000°K  case,  the  radiation  comes  Dredomi- 
nantly  from  the  optically  thin  molecular  bands,  hence  K  is  independent 
cf  s.  Note,  however,  that  the  constant  value  of  K  is  not  that  of  the 
Planck  mean,  due  to  the  atomic  lines  which  are  self-absorbed  and  insigni¬ 
ficant  contributors  to  the  radiation  at  the  conditions  for  which  the 
curves  are  plotted  but  are  not  insignificant  in  an  optically  thin  medium. 
In  fact,  K  does  not  equal  under  any  of  the  conditions  presented, 

which  indicates  that  the  lines  are  never  all  optically  thin  at  these 
conditions. 

An  explanation  for  the  fall  of  K  with  s  for  the  case  when  the 


lines  are  self-absorbed  but  the  continuum  is  optically  thin  (which  is 
the  most  common  situation)  will  now  be  given.  Suppose  one  calculates  the 
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value  of  ft  at  a  particular  value  of  a.  This  value  of  ft  will  always 
correspond  to  optically  thin,  gray  radiation  (see  the  values  of  ft  and  s 
of  Figs.  12b  and  12c),  Thus  for  increases  of  path  length,  the  use  of  a 
constant  value  of  ft  would  predict  a  radiation  intensity  which  grows  lin¬ 
early  with  s  whereas  in  the  real  gas  the  radiation  will  grow  more  slowly 
due  to  self-absorption  of  the  lines.  Thus,  to  match  the  "gray'1  intensity 
to  the  true  intensity,  ft  must  decrease  with  s. 

The  Influence  of  lbppler  Broadening 

Previous  workers  did  not  include  Doppler  broadening  in  their 
calculations  of  line  intensities  since  preliminary  estimates  may  well  lead 
one  to  expect  its  effect  to  be  negligible.  Some  of  our  own  early  calcula¬ 
tions  also  neglected  Doppler  broadening.  However,  it  turns  out  from  our 
more  accurate  calculations  that  Doppler  broadening  has  a  significant  influ 
ence  on  the  amount  of  energy  radiated  at  temperatures  above  approximately 
10,000°K  provided  that  the  density  is  not  too  high.  We  will  therefore  dis 
cuss  the  effect  of  Doppler  broadening  in  some  detail. 

From  the  discussion  of  Part  1,  Section  4,7,  it  can  be  seen  that  the 
influence  of  Doppler  broadening  is  at  its  greatest  when  the  dominant  lines 
are  low-lying  (and  hence  have  significant  Doppler  broadening)  and  when 
these  lines  are  self-absorbed  but  not  strongly  self-absorbed.  (When 
strongly  self-absorbed,  the  intensity  of  a  line  depends  on  the  Stark  width 
and  not  on  the  Doppler  width.) 

We  are  able  to  illustrated  the  quantitative  effect  of  Doppler 
broadtaing  on  the  total  intensity  under  certain  conditions  because  o*  the 
early  calculations  which  were  performed  neglecting  Doppler  broadening. 

The  ratio  of  the  results  from  these  early  calculations  to  the  total  inten¬ 
sity  obtained  from  the  final  calculations  is  plotted  on  Fig.  lu.  Values 
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it 

are  available  at  20,000°K  for  three  densities  and  at  11,000°K  for 

o/o  «  10  It  can  be  seen  from  Tig.  14  that  even,  at  these  moderate  tem- 

o 

peratures,  large  errors  can  be  caused  by  the  neglect  of  Doppler  broadening. 
In  general,  the  error  decreases  with  path  length  as  self-absorption  in¬ 
creases  and  the  controlling  part  of  more  and  more  lines  moves  into  the 
Stark -broadened  wings.  At  sufficiently  lew  optical  depths,  however,  the 
error  increases  with  increasing  s  as  an  increasing  number  of  lines 
become  self  absorbed. 

In  general,  the  influence  of  Doppler  broadening  will  increase  with 
increasing  temperature  and  with  decreasing  density  for  the  following  rea¬ 
sons.  Increasing  temperature  increases  the  Doppler  line  width  and  the 
importance  of  the  Doppler-dominated  lines  associated  with  the  higher  ions. 
Decreasing  the  density  lowers  the  optical  depth  so  that  an  increasing 
number  of  lines  fall  in  the  Doppler-dominated  region  of  their  curves  of 
growth.  At  the  same  time,  decreasing  the  density  increases  the  dominance 
of  the  higher  ions  which  are  predominantly  Doppler  broadened. 


*  -5 

The  values  plotted  for  the  conditions  T  *  20,000°K  and  o/o  »  10 

overestimate  the  correct  values  of  I^/I  since  the  calculation  of  line 

transfer  by  the  third  spectrum  (NIII)  in  the  absence  of  Doppler  broaden¬ 
ing  used  values  of  electron  impact  line  width  which  were  too  large  by  a 
factor  of  10. 
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TABLE  4.1 

Line  Classification  and  Symbols  Used  on  rips.  0  to  11 

First  Spectrum,  NI 


Class  A  Lines 


Croup 

_  c  _  No.  of  Transitions 

Type  of  Transition  _  «  T  j •  •  *  ,, 

Treated  Individually 

Symbol 

1 

2s  -  2p 

6 

IA1 

2 

( 3P)2p  -  ( 3P)3f 

6 

IA2 

3 

31  -  3?.' 

14 

IA3 

4 

31  -  4?.' 

19 

IA4 

5 

4£  -  4?,* 

4 

IA3 

6 

4 9.  -  se* 

6 

IA6 

7 

51-51’,  6*-6fc‘ 

4 

IA7 

Class  B  Lines 

Croup 

_  r  No.  of  Series 

Type  of  Transition  _  ii 

J  Treated  Individually 

Symbol 

1 

(3P)2p-(3P)nf  n  ^  4 

3 

IB1 

2 

(1Df1S)2p-(1D,1S)n«  n  >_  3 

3 

IB2 

f  2s2p*’-2s2p3nfc  , 

L  2s22p3-2s2p3nfc,  "* 

10 

IB3 

4 

3t  -  nl'  n  >  5 

3 

IB4 

5 

4?  -  nl'  n  >  6 

1 

IB5 

Class  C 

contains  all  unmerped  series  with  initial  quantum  numbers  m 

The  symbol  is  IC. 

Second  Spectrum, 

Nil 

Class  A  Lines 

Croup 

,  _  .  .  No.  of  Transitions 

*yoe  of  Transition  _  .  ,  .  ,, 

ireated  Individually 

Symbol 

X 

2s  -  2p 

6 

IIAi 

2 

31  -  3V 

8 

I IA2 

3 

3t  -  41' 

11 

I TA3 

4 

4 1  -  41’ 

2 

IIA4 

Class  B  I ines 

Croup 

.  _  .  .  No. 

of  Series 

Symbol 

»ype  of  iransition  ^  ..  ,, 

3  Treated  Individually 

1 

2p  -  nt  n  ^  3 

3 

TIBI 

2 

31  -  nl'  n>5 

3 

TIB? 

Class  C  contains  all  unmerped  series  with  initial  quantum  numbers  m  >  4. 
The  symbol  is  TIC. 
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Group 


TABLE  4.1  (Continued) 

Third  Snectrum,  .’III 
Class  A  Lines 

Type  of  Transition  of  transitions 

Treated  Indiviuually  bymDoi 

Oc 2 o_2  ~  _ 


2sTp  -  2s2p" 
2s2p2  -  2p3 
^  3s  -  3p  j 
L  4s  -  4p  J 


Class  R  Lines 


Group 

1 


Type  of  Transition  Mo*  of  Scries 

Treated  Individually 

2p  •  nl  n  >  3  1 


IIIA1 

IIIA2 

IIIA3 


Symbol 

III31 


Class  C  contains  all  unmerped  series  with  initial  quantum  nu  ibers  m  >  3. 
Tne  symbol  is  IIIC.  '  — 

Fourth  Spectrum.  NIV 
Class  A  Lines 


Group 

1 


Type  of  Transition  J°'  of  transitions 

Treated  Individually 

'  2s2  -  2s2p  1 

i  2s2n  -  )  4  IVA1 


Group 

1 


i  2s2p  -  2p2  )  4  IVA1 

Class  B  Lines 

Type  of  Transition  _  No*  of  Series  cVP1v,_1 

Treated  Individually  ’ynbol 


2s  -  np  n  >_  3 


IVRl 


Class  c  contains  all  unmerped  series  with  initial  quantum  numbers  m  >  3. 
The  symbol  is  IVC.  — 


Continuum  and  ^seudo-Continuum  Processes 


Process 

Symbol 

NI 

Bound -T ree 

TBF 

NI 

Hiph  Lines 

ILL 

JiT 

iiiph  Geries 

IMS 

nt  r 

Bound-F  ree 

Ilpr 

Nil 

h’iph  Lines 

1 1  ML 

Nil  I 

Bound-F  ree 

T  T  IBr 

NIT  I 

Hiph  Lines 

1 1  TEL 

NT  1 1 

Miph  Series 

TTIHS 

NIV 

Bourd-Free 

IVBr 
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4 

t 


TA^LH  4,1  (Continued) 


Process  Symbol 

Hiph  Lines  TVML 

Hiph  Series  IVHS 

Free-Free  mPF 

Free-Free  2tFF 

Free-Free  3+FF 

Free-Free  4+PF 

Molecular  Bands  2MB 

Molecular  Bands  2+MB 


Note:  ^portion  of  the  continuum  without  a  symbol  shows  the  combined 
e^ect  o*  var minor  contributors  plus  the  influence  of 
seif-absorption  (see  Section  4.3). 
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Tabla  *. J 

Typical  Output  of  Progran  RAO* 


KAL'lATIVE  ENERGY  TRANTER  IN  I  SQTHERNAl  Nl TPOGCN 


SPKIfcS  tC'MUSI  IlfA  At  I'-PlrfATuM*  1SCOC.O  DEGREES  KELVIN,  DENSITY-  O.IOOOE-Of  I  NONO I  MENS  ION  AL  I 

Co.  RE  SPONGING  P«E  SSUftE  »  0.1S931  00  ATNOSPMERES 


values  m 

pahticle  density  Oiviofr 

BY  IOSCHPI0T 

NUNBE*  l-f. 

fcRTtElNi 

K  /  • 

Nf 

N 

N* 

NM 

N*M 

N*M* 

NE 

c.iEU.’r-r» 

t  :.v/iue-c> 

0. I6f RNE-C? 

0. ?  JT  TAE-07 

O.IOfAlE-lB 

O.ASiSNE- 

0.  I6289f-0f 

PART ITICN 

E UNCTIONS 

Of* 

CN  2 

CN 

UN* 

CNM 

as»»» 

ON*  ♦♦ ♦ 

l.  1,’cVM  ct. 

"J 

J 

n 

c 

r 

o.iofsn  cf 

o.oonvu  oi 

0. 101 AfE  01 

O.fOOfbE 

01 

GAS  CCAUlTIONS 

Tl  -PCRAURf  -  ISCCC.O  DEGREES  KELVIN 

DENSITY,  O.lOOOf-Of  IVONOI  KENS  IfJNAL  I 
RATH  IfNCTE*  0.1000E  Of  C*>. 


0»M»AI  AfSUlTS 


U:ial  INI(NSITy> 


ILL«  ►«C"  '  '|V  f*  H(«MSS  O.'.v,*  L I  » 

inunsi  tv  n  e  ic  i  im 
ISttHWT*  CUE  I '  Ci'M  I  Vjl!»« 
LIM  1-TtASIIY  All"  -At-CAIUM  CONTINUE* 
T'lVl'TIC  AY 

i  im  i»uvsiit  •.•iC*r»«rcAO  continui*  ■ 

UllYAlrM  l«*Y-r.AS  A-t\^at»m\  CCt  EE  IC  It  N  T  ■ 


O.S/fVE-OA 
0.139TE  H 
O.lSTlf  -0  1 
.  .*  J90E  0  A 
J. ANA  II  ■  '»  C 
C.AJAVE  OC 


fRG  (V  C»-f  STfRAO-l 
ERG  SIC-I  CN-f  STEtAO-l 

»  BL ACA-BO0Y  intensity 
yams  C*W 

»  IOTAI  INTENSITY 
i  TOTAL  INTENSITY 


vi.9fO»E  OC 
O.INAIE-O/  C»-5 
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CCNTINUUK  INTENSITIES  RELATIVE  TO  THE  TOTAL  INTENSITY 

ALL  CONTRIBUTORS  0.535SE  00 
N  f  BOUND-FREE  0.3844E  00 

NI  HIGH  LINES  0.I205E  00 

Nil  BOUNO-EREE  0.0 

MI  HIGH  LINES  0.0 

Nil  I  BOUND-FREE  C.O 
N 1 1 1  HIGH  LINES  C.O 
NIV  BOUND-FREE  0.0 
NIV  HIGH  LINES  C.O 

K+  FREE-FREE  0.2651F-01 
NI  HIGH  SERIES  0.0 

N*»  FREE-FREE  0.0 

Nil  HIGH  SERIES  0.0 

N^**  FREE-FREE  C.O 
Nil  I  HIGH  SERIES  0.0 
*++*+  FREE-FREE  C.O 
NIV  HIGH  SERIES  0.0 
N2  MOLECULAR  BANOS  0.0 
N2*  MOLECULAR  BANOS  0.0 
ALL  NI  CONTRIBUTORS  C.50*9E  00 
ALL  Nil  CONTRIBUTORS  C.2651E-01 
ALL  Mil  CONTRIBUTORS  0.0 
ALL  NIV  CONTRIBUTORS  C.O 

SUM  CF  INL I VI  DUALLY  CALCULATED  CONTRIBUTORS  0.5315E  00 


NOTE:  VALUES  FOR  INDIVIDUAL  CONTRIBUTORS 
ARE  ACCURATF  TO  APPROX  I MATELY  ♦OR  -  0.I1E-0! 
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LIKE  INTENSI T i FS 


HE l AT l Vt  TO  THE  TCTAl  INTENSI T 

W 

RELATIVE  TO  INTENSITIES 

IN  THE  PRESENCE  UF 

IN 

THE  ABSENCE  OE 

A  MCKCPCUNC  CCNTINUUH 

A  BACKGROUND  CONTINUUM 

f  IRST 

SPECTRUM 

ALL 

LINES  RELATIVE  TC  TOTAL 

INTENSITY 

0.45IOE  CO 

CLASS 

A  LINES 

I  0.9072E-02 

1 

0.9269E  00 

2  0.  MS5C-CI 

2 

0.9798E  00 

3  C.Ss>62E-Cl 

3 

0.995BE  00 

4  C.12C1F-0I 

6 

0.9907E  00 

S  O.R6b7t-03 

5 

0.9626E  00 

6  C.520RE-C1 

6 

0.9420E  00 

7  0.6325E-04 

7 

0. 7 BO°E  Ou 

CLASS 

5  LINES 

GROUP  INTENSITY 

GROUP 

INTENSITY 

1  C.l«5SE  CC 

1 

0.9174E  00 

i  or, 

2 

0.8708E  00 

3  C.VOCF-02 

3 

0.9102E  00 

A  C.O 

4 

0.0 

5  C.O 

5 

0.0 

CLASS 

C  LINES 

IMENSI  TY 

INTENSITY 

C.R*74t-C3 

0.8069E 

00 

SECCNIJ 

SPECTRUM 

*LL 

L  I*FS  RELATIVE  TO  TOTAL 

INTENSI TY 

0. 13O.E-01 

CLASS 

A  LINES 

i  c.msE-ci 

l 

0.906 7E  00 

2  C.3SO3E-03 

2 

0.9972E  CO 

3  0.  77B4F  -0*. 

3 

0.99A0F  00 

*.  0.12  ISP-05 

4 

0.9909E  0? 
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MISCELLANEOUS  LINE  INFORMATION 


THIRO  SPECTRUM 
FIRST  MERGED  SERIFS  9 

FIRST  MERGED  LINE 
Of  NCNHYQRCGEN I C  SERIES 
GROUP  SERIES  LINE 
I  1  6 

FIRST  MERGEC  UNE 
OF  HYDRGGENIC  SFRIES 
M  N 

3  9 

4  9 

5  9 

6  9 

7  9 

8  9 

FOURTH  SPECTRUM 
FIRST  MERGED  SERIFS  9 
FIRST  MERGEC  LINE 
OF  NONHYDRCGENIC  SERIFS 
GROUP  SERIES  LINE 
l  I  6 

FIRST  MERGED  LINE 
OF  HYDROGEN  I C  SERIFS 
M  N 

3  9 

4  9 

5  9 

6  9 

7  q 

F  9 


SECOND  SPECTRUM 
FIRST  MERGED  SSRIES  9 
FIRST  MERGEC  LINE 
CF  NONHYOROGE'NIC  SERIFS 
GROUP  SERIES  LINE 
1  i  II 

l  2  10 

1  3  6 

2  I  5 

2  2  5 

2  3  5 

FIRST  MERGED  LINE 
OF  HYCROGENIC  SERIES 
M  N 

4  9 

5  9 

6  9 

7  9 

8  9 
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TABLE  4.3 


A  comparison  of  the  Calculated  ^lanck  Mean  Absorption 

7 

Coefficient  with  the  Results  of  Armstrong  et  al. 


Planck  Mean  Absorption  Coefficient  cm 


kT  eV 

o/o 

Present 

Calculations 

Ref 

7  1 
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0.98  «  io_1  : 

1.37 

1.74 

"  j 
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10-J 

1.35  x 

io-1 

1.71  x 

io-1  I 

io-3  ; 

1.16  x 

io-2 

1.54  x 

io'2  ' 

10-“ 
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0,78  x 

.io-3 
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10-5 
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APPENDIX  T  A  Derivation  of  the  Directionally  Approximate  Flux  Equations 

Many  derivations  of  directionally  approximate  equations  of  radiative 

30  31 

transfer  may  be  found  in  the  literature  *  .  We  are  interested  here  in 

one-dimensional,  half-range  flux  equations.  Tt  does  not  appear  to  have  been 
pointed  out  previously  that  such  equations  can  be  derived  with  two  free 
constants  by  applying  the  well-known  substitute  kernel  approximation  to 
the  half-range  fluxes.  This  will  now  be  demonstrated. 

The  half-range  fluxes  q*  and  q»~  at  optical  thickness  are 
defined  as  follows 


/•l 

yV  - 2’  V  d" 

*  n 


(1.1) 


and 


qv  =  2n 


r  o 


-1 


IyM  du. 


(1.2) 


A  formal  solution  of  the  equation  of  radiative  transfer  for  the  half-range 
fluxes  in  a  one-dimensional  slab  with  isotropic  wall  conditions  yields 


4-  fV 

%  =  2<C(o)W  * 2ir  Vy)  Vnvy)dy 

*  n 


(1. 3a) 


and 


V’t)E3(Vnv> 


2*  By(y)  E2(y-nv)dy 

*  n 


(I. 3b) 


where  E  (x)  =  e'"x^1  un”2du  (see  Ref.  32),  Making  the  substitute  kernel 

n  o 

%  — bx  ^ 

approximations  E2(x)  *  ae  ano  E^(x)  *  co  and  dif ferentiating 
yields  the  directionally  approximate  equations 
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and 


s  2iraB 

v 


(1. 4a) 


dq  _ 

-  *  2iraB  +-bq".  (1. 4b) 

dn  v  v 

v 


As  is  well  known,  Eqs.  (I. 4a)  and  (I. 4b)  can  be  combined  to  produce  a 
second  order  equation  in  the  full-range  flux. 

Eqs.  (I. 4a)  and  (I. 4b)  have  two  free  constants  a  and  b.  There  are 
a  number  of  ways  of  choosing  these  constants  but  we  will  mention  only  one; 
the  constants  may  be  chosen  to  give  the  correct  expressions  for  half-range 
flux  in  the  two  limits  of  optical  depth  approaching  zero  and  infinity. 
Matching  Eqs.  (I. 4a)  and  (I. 4b)  to  the  exact  expressions  in  these  two 
limits  yields  a  *  1  and  b  *  2.  (Note  that  these  values  do  not  satisfy 
the  optically  thick  limit  of  the  full  range  flux  equation.) 
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APPENDIX  II  Fortran  Program  RADN 

This  appendix  gives  a  brief  description  of  the  program  RADN  which 
calculates  the  specific  intensity  at  a  point  in  uniform  nitrogen  gas.  The 
main  program  (containing  630  statements)  and  16  subprograms  (containing  a 
total  of  920  statements)  which  comprise  RADN  are  discussed  in  turn.  Each 
was  written  in  FORTRAN  IV  (E-level  subset)  for  the  IBM  360  computer  (Model 
50)  at  Brown  University. 

Program  RADN  takes  from  3  to  10  minutes  to  calculate  the  results  for 
a  single  point  (given  temperature,  density  and  path  length),  depending  on 
the  conditions  and  the  number  of  individual  continuum  contributions 
investigated. 

Main  Program  RADN 

The  main  program  reads  data  for  the  composition  calculation,  data  for 
Class  B  lines  of  HI,  Nil,  NIII  and  NIV,  data  for  Class  A  lines  of  NI,  Nil, 
NIII  and  NIV,  data  for  bound-free  transitions  of  NT,  KII,  Nil I  and  NIV, 
one  value  of  temperature,  data  for  molecular  bands  (if  the  temperature  is 
5,000°K,  8,000°K  or  11,000°K  only),  up  to  ten  values  of  density,  up  to  ten 
values  of  path-length  and,  finally  a  set  of  parameters  for  each  density 
which  enable  a  selection  of  the  individual  continuum  contributors  to  be 
omitted. 

The  main  program  then  calls  and  organizes  the  subprograms  to  calcu¬ 
late  line  and  continuum  intensities.  Finally,  it  processes  and  writes  the 
results. 

Subprograms  COMPN,  SAHA  and  SUHfi 

These  three  subprograms  calculate  the  composition  according  to  the 
method  described  in  Part  1,  Chapt.  2. 
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Subprogram  MULT 

Subprogram  MULT  calculates  the  intensity  due  to  a  group  of  closely 
spaced  lines  and  indicates  whether  the  group  is  optically  thick  or  optically 
thin.  Overlapping  of  the  lines  is  allowed  for  as  described  in  Part  1,  Sec¬ 
tions  3.4  and  4.7.  The  line  wings  have  either  dispersion  or  quasi-static 
form  and,  where  instructed  to,  the  rubprogram  allows  for  a  Doppler 
broadened  core. 

Subprogram  SOLE 

Subprogram  SOLE  calculates  the  intensity  due  to  each  of  several  groups 
of  (Class  A)  lines.  The  bulk  of  the  data  is  given  in  a  two-dimensional 
array.  For  each  transition,  the  subprogram  will  calculate  the  occupation 
number,  the  Stark  line  width  and,  where  necessary,  the  resonance  line 
width.  These  and  other  necessary  values  are  then  fed  to  subprogram  MULT 
which  returns  the  corresponding  value  of  intensity  (MULT  may  allow  for  a 
Doppler  core  if  the  line  data  instructs  it  to).  SOLE  then  calls  the 
continuum  subprograms  and  finds  the  total  continuum  absorption  coefficient, 

K  at  the  line  center  and  calculates  an  intensity  reduced  by  the  factor 

vc 

exp(-K  s).  Reduced  and  unreduced  intensities  are  summed  for  each  group 
o*c 

separately  and  the  sums  passed  back  *o  the  main  program. 

Subprogram  SRI ESI 

Subprogram  SFIES1  finds  the  position  of  the  pseudo-bound-free 
threshold  for  a  number  of  non-hydrogenic  series.  Each  position  is  stored 
according  to  a  given  parameter  which  enables  the  position  to  be  related  to 
the  correct  photoionization  transition. 

Data  for  the  non-hydrogenic  (Class  B)  lines  is  given  in  a  three- 
dimensional  array,  for  each  transition,  starting  with  the  lowest  lying. 
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the  subprogram  calculates  the  occupation  number  and  Stark  line  width.  These 
and  other  values  are  fed  to  subprogram  MULT  which  returns  the  corresponding 
value  of  intensity  and  indicates  the  optical  thickness  of  the  line.  For 
the  line  to  be  considered  the  start  of  the  pseudo-continuum,  it  must  be 
higher  up  the  series  than  a  specified  minimum  and  be  either  optically  thin 
or  the  first  of  three  lines  which  are  merged  according  to  their  effective 
widths  (see  Part  2t  Section  3.3).  (It  is  necessary  to  check  three  lines 
to  avoid  the  case  where  an  error  in  the  data  results  in  two  lines  being 
closer  together  than  they  should  be.)  Transitions  are  only  considered  for 
merging  if  the  entire  group  of  closely  spaced  lines  is  itself  merged  and 
behaving  like  a  superline  (see  Part  1,  Section  3.4). 

The  intensities  of  all  lines  not  counted  in  the  pseudo-continuum  are 
added  in  groups  and  the  sums  passed  to  the  main  program.  As  in  subprogram 
SOLE,  the  lines  are  stored  both  unreduced  and  reduced  by  a  background  con¬ 
tinuum.  In  this  case,  however,  the  high  lines  are  not  included  in  the 
continuum  (because  the  extent  of  the  high  lines  is  found  within  the  sub¬ 
program  which  would  havs  to  be. called  twice  for  it  to  be  possible  to  include 
high  lines  in  the  background  continuum). 

Subprogram  SRIES2 

Subprogram  SRIES2  performs  operations  nimilar  to  those  executed  by 
SRIES1  but  it  treats  hydrogenic  (Class  C)  series.  Since  the  hydrogenic 
data  can  be  calculated  within  the  subprogram,  only  overall  quantities  such 
as  ionization  potentials  and  core  statistical  weights  need  be  supplied. 

Subprogram  KEFHE 

Subprogram  KEF^E  finds  the  first  completely  merged  or  thin  series 
above  a  specified  minimum  (see  Part  2,  Section3.5).  It  treats  the 


I 


89 


series  as  hydrogenic  and  quasi-statically  broadened. 

Subprogram  C1BF 

Subprogram  C1BF  is  the  first  of  6  subprograms  for  finding  the 
absorption  coefficient  due  to  continuum  processes.  C1BF  finds  the  absorp¬ 
tion  coefficient  due  to  NT  true  bound-free  transitions.  It  requires  some 
cross  sections  in  tabular  form  (we  use  the  Burgess-Seaton  calculations  of 
Sherman  and  Kulander^,  see  Section  3.2)  at  the  threshold  frequencies  which 
must  also  be  given.  Intermediate  values  are  found  by  linear  interpolation. 
Occupation  numbers  are  calculated  in  the  main  program  (from  given  initial 
energies)  and  are  supplied  as  data  to  the  subprogram.  For  higher  initial 
quantum  numbers,  hydrogenic  cross  sections,  occupation  numbers  and 
threshold  frequencies  are  calculated. 

Subprogram  C2HI 

Subprogram  C2HI  calculates  the  absorption  coefficient  due  to  the 
merged  or  thin  high  lines  of  NI.  It  requires  the  same  data  as  C1BF  plus 
the  pseudo-threshold  frequencies  found  by  SRIES1.  For  each  initial  quan¬ 
tum  number  for  which  the  given  frequency  lies  between  a  pseudo-threshold 
and  a  true  threshold,  a  value  of  absorption  coefficient  is  added  to  a  sum. 
If  the  initial  state  is  low-lying,  the  cross  section  is  taken  equal  to  the 
threshold  bound-free  value  (see  Section  2.4).  For  higher  values  it  is 
calculated  from  the  hydrogenic  bound-free  formula. 

Subprogram  C3BF 

Subprogram  C3BF  calculates  the  absorption  coefficient  due  to  ionic 

bound-free  transitions.  The  treatment  is  similar  to  that  of  C1BF  except 

that  the  low-lying  cross  sections  are  computed  from  a  formula  due  to 
22 

Seaton  (see  Section  3.2). 
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Subprogram  C4HI 

Subprogram  C4HI  calculates  the  absorption  coefficient  due  to  ionic  high 
lines*  The  treatment  is  similar  to  that  of  C2HI  except  that  Seaton's  form¬ 
ula  is  used  for  the  low-lying  cross  sections. 

Subprogram  C5FT 

Subprogram  C5FF  calculates  the  absorption  coefficient  due  to  free-free 
transitions  and  high  series  (see  Section  3.3).  Four  parameters  enable  any 
combination  of  these  contributors  to  be  included. 

Subprogram  C6MB 

Subprogram  C6MB  calculates  the  absorption  coefficient  due  to 
molecular  bands.  This  facility  is  available  at  5f000°K,  8„Q00°K  and 
llf000°K  only.  Details  of  the  bands  and  their  treatment  may  be  found  in 
Section  3.4.  It  is  possible  to  omit  either  all  the  bands  of  or  all 

the  bands  of  if  desired. 

Subprogram  SIMP 

Subprogram  SIMP  performs  a  Simpson's  rule  integration  on  the  quantity 
B  [l-exp(-K  s)]  from  hv  to  hv  .  It  proceeds  as  follows.  Starting 

V  V  |C  X  a 

with  N  intervals,  it  keeps  doubling  the  number  of  intervals  until  the 
fractional  error  in  two  successive  integrations  is  less  than  a  given  quan¬ 
tity.  If  this  match  is  not  achieved  in  6  attempts,  the  subprogram  accepts 
the  last  value  but  writes  a  warning  message  stating  the  oustanding  error, 
the  maximum  error  at  which  a  match  would  have  been  accepted  and  the  value 
of  a  parameter  which  enables  one  to  identify  the  continuum  process  being 
investigated  when  the  unmatched  integration  occurred.  The  values  of  K 

v,c 

are  computed  in  this  subprogram  by  means  of  the  continuum  subprograms 


(CIBr  to  C6MB) . 
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Subprogram  CONT 

Subprogram  CONT  calculates  the  intensity  due  to  continuum  processes. 

From  0  to  hv.  and  from  hv  to  •.  the  intensity  is  calculated  as 

described  in  Section  4,2.  From  hv^  =  0.35kT  to  hv^  =  8,9kT,  the  inten 

sity  is  calculated  by  SIMP  with  N  =  50  to  a  match  of  1%;  this  result  we 

call  FINE.  From  hv.  to  hv,  and  from  hv„  to  hv  .  the  photon 

lr  1  2  uv*  y 

energy  is  divided  into  a  number  of  intervals  hv.  to  hv.  ,  where 

3  j+1 

v^+i  *  100  *  Vj.  Each  interval  is  calculated  by  SIMP  with  N  =  10  to  a 
match  of  10%.  If  the  integrated  intensity  of  any  interval  exceeds  5%  of 
FINE  then  the  integration  over  that  interval  is  recomputed  with  N  »  50 
to  a  match  of  1%.  In  this  way,  an  overall  accuracy  of  approximately  1% 
is  maintained. 
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FIG.  I  SCHEMATIC  DIAGRAM  SHOWING  SOME  OF  THE 
ENERGY  LEVELS  OF  ATOMIC  NITROGEN 
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IN  ADDITION  TO  THOSE  SHOWN, THE  FOLLOWING 
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FIG.  2  SCHEMATIC  DiAGRAM  SHOWING  THE  MOST  IMPORTANT  ENERGY  LEVELS  OF 
ATOMIC  NITROGEN 
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IN  ADDITION  TO  THOSE  SHOWN.  THE  FOLLOWING 
CONFIGURATIONS  OCCUR 


CONFIGURATION 

IONISATION  EDGE  (ev ) 

2pVs)ni 

52.7 

2ps{lD)ni 

54  8 

2p*(*P)ni 

58  2 

FIG.  3  SCHEMATIC  DIAGRAM  SHOWING  THE  MOST  IMPORTANT  ENERGY  LEVELS  OF 
SINGLY  ION L  ED  NITROGEN 


98 


IN  ADDITION  TO  THOSE  SHOW  N  ,THE  FOLLOWING 
CONFIGURATIONS  OCCUR  : 
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FIG.  4  SCHEMATIC  DIAGRAM  SHOWING  THE  MOST  IMPORTANT 
ENERGY  LEVELS  OF  DOUBLY  IONISED  NITROGEN 
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Notes  relevant  to  Pig.  5: 

Each  horizontal  line  shows  the  range  of  frequencies  of  a  class  of  spectral 
lines. 

For  each  class  of  spectral  lines,  the  type  of  initial  state  is  shown  to  the 
left  of  the  plotted  line. 

Where  parentage  splitting  can  occur,  the  core  term  is  shown  in  parentheses 
to  the  right  of  the  plotted  line. 

Both  lower  and  upper  states  are  shown  for  same-shell  transitions. 

Within  each  spectrum,  the  plotted  lines  are  arranged  in  vertical  order 
approximately  according  to  the  energy  of  the  initial  state.  Apart  from 
this  arrangement,  the  vertical  distribution  has  no  significance. 


FIG.  6a  RATIO  OF  FREQUENCY  INTEGRATED  INTENSITY 
TO  BLACK-BODY  INTENSITY  AS  A  FUNCTION  OF 
PATH  LENGTH 
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FIG. 7  RATIO  OF  FREQUENCY  INTEGRATED  INTENSITY 
TO  BLACK -BODY  INTENSITY  DIVIDED  BY  DENSITY 
AS  A  FUNCTION  OF  TEMPERATURE  FOR 
PATH  LENGTH  *  10  cm 
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A  key  to  the  symbols  of  Fig9.  9  to  11 
may  be  found  in  Table  4.1 
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FIG.  9a  DISTRIBUTION  OF  RADIATIVE  INTENSITY 
OVER  INDIVIDUAL  CONTRIBUTORS  AS  A 
FUNCTION  OF  PATH  LENGTH  FOR 
TEMPERATURE,  T=5,000°K 
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FIG.  9b  T  =  II , 000 °K 
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FIG.  9c  T *  20,000  °K 
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FIG.  9e  T»  35 , 000 °K 
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FIG.  10  DISTRIBUTION  OF  RADIATIVE  INTENSITY 
OVER  INDIVIDUAL  CONTRIBUTORS  AS  A 
FUNCTION  OF  TEMPERATURE  FOR 
PATH  LENGTH,  s  =  IOcm 
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T  s  11,000  °K 


FIG. II  DISTRIBUTION  OF  RADIATIVE  INTENSITY 
OVER  INDIVIDUAL  CONTRIBUTORS  AS  A 
FUNCTION  OF  DENSITY  FOR 
PATH  LENGTH  .SMOClTl 
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VALUES  OF  PLANCK  MEAN  ABSORPTION  COEFFICIENT 
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VALUES  OF  PLANCK  MEAN  ABSORPTION  COEFFICIENT 

TtfK  5,000  11,000  15,000  20,000  30,000  35,000 

KpCfTT1  2.28x10"*  3.62XI0'5  l.88x|Q-4  5.65xl0"4  1.30x10-*  I.SSxlO"3 
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FIG.  13  DISTRIBUTION  OF  PARTICLE  NUMBER 

DENSITIES  OF  RADIATING  SPECIES  OVER 
INDIVIDUAL  CONTRIBUTORS  AS  A 
FUNCTION  OF  TEMPERATURE. 

NOTE:THE  PARTICLE  NUMBER  DENSITY  OF 
N j"  IS  NEVER  SUFFICIENTLY  HIGH  TO  APPEAR 
ON  THIS  REPRESENTATION 
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FIG.  14  THE  INFLUENCE  OF  DOPPLER 
BROADENING  ON  THE  TOTAL 
INTENSITY  FROM  NITROGEN  GAS. 
IND  IS  THE  TOTAL  INTENSITY 

WHEN  DOPPLER  BROADENING  OF 
THE  LINES  IS  OMITTED. 

X  IS  THE  TOTAL  INI  ENSI  TY  WHEN 
DOPPLER  BROADENING  OF  THE 
LINES  IS  INCLUDED. 
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